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Abstract 

A workable  open  loop  guidance  scheme  for  orbital  transfer 
maneuvers  is  developed  for  a two  stage  solid-rocket  vehicle 
which  has  no  thrust  termination  capability.  The  scheme  ef- 
fectively manages  any  excess  energy  by  matching  a non-Hohmann 
transfer  trajectory  to  the  fixed  energy  (aV)  capabilities  of 
the  vehicle.  The  entire  burden  of  effecting  the  transfer  is 
put  on  prelaunch  targeting,  so  that  during  the  bums  the  thrust 
can  be  directed  along  a precomputed  direction  using  constant 
attitude  maneuvers  only. 

A computer  program  has  been  developed  which  employs  a 
nonlinear  equation  solving  routine  to  accomplish  exact  tar- 
geting for  the  finite-thrust  transfer  maneuver.  The  transfer 
trajectory  is  characterized  by  six  control  parameters  (the 
outputs  of  targeting),  and  the  final  orbit  is  defined  by  a 
set  of  "hit  conditions".  The  values  of  the  control  parameters 
which  drive  the  vehicle  state  vector  to  satisfy  the  hit  con- 
ditions become  the  guidance  system  target  parameters. 

In  addition,  an  error  analysis  is  performed  on  the  scheme 
throughout  the  range  of  possible  trajectories  which  exist  for 
excess  energy  missions.  These  trajectories  are  then  compared 
on  the  basis  of  optimality,  such  as  minimum  insertion  errors 
and  transfer  time.  Results  are  presented  for  geosynchronous 
and  subsynchronous  transfers  between  circular  orbits. 
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AN  ENERGY  MANAGEMENT  GUIDANCE  SCHEME  APPLICABLE 


I 


TO  THE  INTERIM  UPPER  STAGE 


I . Introduction 

Background 

The  United  States  space  program  is  currently  focused  on 
the  development  and  implementation  of  the  Space  Transportation 
System  (STS),  better  known  as  the  "Space  Shuttle".  This  is  to 
be  a system  that  will  serve  the  routine  operational  space  re- 
quirements in  the  1980  decade  and  beyond. 

The  major  component  of  the  space  shuttle  is  the  "Orbiter" 
vehicle,  which  somewhat  resembles  a cargo-type  aircraft.  The 
orbiter  will  be  boosted  into  low  earth  orbit,  and  after  com- 
pleting its  mission,  well  re-enter  the  atmosphere  and  glide 
to  a landing  much  like  a conventional  powered  aircraft.  Due 
to  its  size  and  weight,  the  orbiter  vehicle  can  be  placed  in 
only  a relatively  low  earth  orbit,  but  will  carry  extra  pro- 
pulsive stages  in  its  cargo  bay  to  complete  higher  energy  mis- 
sions. 

While  in  orbit,  these  extra  propulsive  stages  will  be 
placed  outside  the  orbiter  where  they  can  be  launched  to  de- 
liver a satellite  or  other  payload  to  a higher  orbit.  The  low 
orbit  is  usually  referred  to  as  the  "parking"  orbit,  and  the 
higher  orbit  as  the  "mission"  orbit.  The  extra  propulsive 
stages,  which  complete  the  orbital  transfer,  are  usually  re- 
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ferred  to  jointly  as  the  "upper  stage"  vehicle. 

The  United  States  Air  Force  Space  and  Missile  Systems 
Organization  (SAMSO)  has  been  tasked  with  the  development  of 
an  Interim  Upper  Stage  (IUS)  vehicle,  which  will  be  used  un- 
til a fuller  capability  vehicle  can  be  designed  and  produced. 
The  "Burner  II"  space  booster,  made  by  the  Boeing  Aerospace 
Company,  has  recently  been  selected  by  SAMSO  for  modification 
and  use  as  the  IUS  vehicle. 

The  baseline  Burner  II  vehicle  consists  of  two  stages, 
where  the  first  stage  burn  is  used  to  place  the  IUS  into  an 
elliptical  transfer  orbit,  and  the  second  stage  bum  is  used 
to  insert  the  IUS  into  the  mission  orbit.  The  Burner  II  uses 
a propulsion  system  consisting  of  two  solid  propellant  rocket 
motors  which  have  no  thrust  termination  capability. 

Since  the  engines  cannot  be  shutdown  prior  to  depletion 
of  all  the  propellant,  and  off-loading  of  solid  propellant 
to  tailor  each  mission  to  its  minimum  energy  requirements  is 
impractical,  the  solid  rocket  motors  will  usually  produce  more 
energy  (velocity  change  capability)  than  is  necessary  to  com- 
plete the  orbital  transfer.  Therefore,  any  guidance  scenario 
used  to  complete  the  transfer  must  involve  some  method  of  de- 
pleting (or  somehow  utilizing)  the  excess  energy.  This  process 
is  termed  "energy  management". 

In  summary,  the  overall  guidance  and  navigation  problem  of 
the  IUS  is  primarily  to  complete  the  orbital  transfer!  but  this 
is  complicated  by  the  requirement  of  managing  any  excess  energy 
in  the  process. 


Guidance.  Targeting  and  Energy  Management 

The  meanings  of  these  three  terms  should  first  be  made 
clear.  Targeting  consists  of  computations  done  prior  to  launch, 
usually  in  a ground  based  computer,  to  supply  mission  dependent 
parameter  values  to  the  on-board  flight  program,  where  they  are 
stored  for  use  during  the  maneuver.  Guidance . in  the  strictest 
sense,  usually  means  on-board  computations  carried  out  in  closed 
loop  fashion  during  the  actual  thrusting  portions  of  the  maneu- 
ver to  provide  steering  commands  for  the  vehicle  propulsion  sys- 
tem. When  the  term  targeting  is  used  in  this  study,  its  mean- 
ing will  adhere  to  the  above  definition  in  a strict  sense.  The 
term  guidance,  however,  will  often  be  used  more  loosely,  and 
will  tend  to  infer  any  and  all  computational  processes  neces- 
sary to  effect  the  orbital  transfer  maneuver. 

The  concept  of  energy  management  was  briefly  introduced 
in  the  last  section.  There  are  many  possibilities  available 
for  handling  the  excess  fuel,  both  in  the  premission  target- 
ing and/or  during  on-board  guidance  phases.  The  fuel  deple- 
tion problem  is  relatively  new,  but  some  work  has  been  done 
on  this  concept  recently.  Several  good  examples  are  early 
proposals  made  by  the  Boeing  Company  for  the  Burner  II  (Ref  2), 
and  by  the  Charles  Stark  Draper  Laboratory  for  the  Navy's 
Trident  Missile  (Ref  3)* 

Both  of  these  proposals  utilized  maneuvering  during  bums 
to  deplete  the  excess  propellant.  The  Boeing  proposal  used  an 
attitude  modulation  technique  that  would  rotate  the  thrust  vec- 
tor to  equal  angles  each  side  of  a nominal  thrust  direction. 

The  Trident  scheme  rotated  the  thrust  vector  ^through  an  arc 
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where  the  arc  length  is  equal  to  the  velocity  change  capabil- 
ity of  the  motor,  and  the  chord  length  is  equal  to  the  veloc- 
ity change  necessary.  Both  ideas  are  shown  in  Figure  1-1. 

The  underlying  objective  of  guidance  (loose  connotation 
again  emphasized)  is  generally  to  place  a vehicle  on  some  form 
of  free  fall  trajectory  which  satisfies  given  specifications 
or  constraints.  Here  the  free  fall  trajectory  to  be  satisfied 
is  the  specified  mission  orbit.  To  completely  specify  an  orbit, 
a maximum  of  six  parameters  (constraints)  must  be  satisfied. 
These  parameters,  which  completely  describe  the  orbit,  are  nor- 
mally chosen  to  be  the  classical  orbital  elements  (Ref  1:58). 


Attitude 

Angle 

/ “A 

(nominal) 

/ \ jime 

Xnr 

Burner  II 

Vg  direction 

Trident 

Figure  1-1.  Prior  Energy  Management  Proposals 
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In  most  guidance  schemes  (with  thrust  cut  off  capability) 
the  idea  is  normally  to  effect  thrusting  in  the  direction  of 
the  desired  velocity  vector.  This  vector  is  called  the  ve- 
locity-to-be-gained  (Vg),  and  is  defined  at  any  instant  of  time 
to  be  the  vector  difference  between  the  velocity  required  (at 
that  instant  to  satisfy  final  constraints),  and  the  actual  ve- 
hicle velocity.  The  usual  method,  then,  is  to  thrust  in  the 
direction  of  the  Vg  vector  in  order  to  drive  it  to  zero  as 
soon  as  possible  (i.e.,  with  the  minimum  expenditure  of  fuel). 
The  instant  Vg  goes  to  zero  the  engine  is  shut  down. 

The  main  point  is  that  thrust  cut  off  capability,  which 
controls  the  magnitude  of  the  velocity  change,  is  an  important 
control  variable  usually  available  for  velocity-to-be-gained 
guidance.  With  thrust  cut  off  capability  there  are  eight  de- 
grees-of-freedom  available  to  satisfy  mission  constraints;  as- 
suming that  the  bums  are  constant  attitude  and  directed  at 
making  up  Vg.  These  eight  degrees-of-freedom  (mission  varia- 
bles) are  the  two  ignition  times,  and  the  three  components  of 
the  velocity  change  (AV)  vectors  for  each  stage.  The  three 
components  of  AV  can  be  thought  of  as  a magnitude,  a pitch 
angle  and  a yaw  angle. 

When  there  is  no  thrust  cut  off  capability  two  degrees- 
of-freedom  are  lost,  since  the  magnitudes  of  each  AV  vector 
are  now  fixed  by  the  total  amount  of  propellant  on-board  each 
stage.  Thus,  for  the  IUS  orbital  transfer  problem,  six  mis- 
sion variables  are  available  (two  ignition  times  and  two  thrust 
direction  angles  for  each  bum),  which  are  sufficient  to  sat- 
isfy a total  of  six  mission  constraints i where  most  or  all  of 
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these  constraints  could  be  the  orbital  parameters  which  specify 
the  desired  mission  orbit. 

In  comparing  energy  management  guidance  with  convention- 
al velocity-to-be-gained  guidance,  the  former  becomes  more 
constrained  in  that  since  the  engine  cannot  be  cut  off,  it  is 
necessary  that  the  fuel  is  somehow  depleted  at  the  exact  time 
the  Vg  goes  to  zero  (i.e.,  the  required  velocity  is  attained). 

An  important  final  point  to  be  made  concerning  any  guidance 
scheme  is  that  the  scheme  must  be  able  to  satisfy  accuracy  re- 
quirements. That  is,  the  overall  navigation,  guidance  and  con- 
trol system  of  the  IUS,  in  whatever  form  it  takes,  must  be  such 
that  errors  in  the  desired  position  and  velocity  vectors  after 
insertion  into  the  mission  orbit,  are  acceptably  small. 

So  in  an  attempt  to  design  any  guidance  scheme,  it  is  of 
overall  importance  to  know  how  errors  propagate  through  the 
scheme  (and  the  resulting  maneuver)  to  produce  errors  after 
insertion  into  the  mission  orbit.  This  is  a practical  measure 
of  the  worth  of  the  scheme,  and  one  of  the  deciding  factors  in 
consideration  of  that  scheme  for  actual  implementation. 

Outline  of  the  Problem 

The  IUS  vehicle  is  to  be  expendable,  so  cost  effective- 
ness is  a most  important  consideration.  Thus,  low  cost  sys- 
tem (hardware  and  software)  requirements  for  guidance  and  con- 
trol would  be  desirable. 

Simplicity  and  reliability  are  also  important.  Simplic- 
ity is  particularly  desirable  to  aid  in  understanding,  and  to 
minimize  both  the  hardware  and  software  requirements.  These 
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points  tend  to  lay  the  ground  rules  and  guidelines  for  this 
study. 

The  goal  was  to  devise  an  energy  management  guidance 
scheme  applicable  to  the  space  shuttle  IUS,  which  could  be 
relatively  simple,  practical  and  cost  effective.  The  scheme 
was  to  have  the  dual  capability  of  completing  the  transfer 
using  the  fixed  velocity  capabilities  of  a two  stage  IUS,  and 
using  a transfer  trajectory  which  would  be  "optimal"  in  some 
sense,  depending  on  the  objectives  of  the  particular  mission. 
An  essentially  "open  loop"  scheme  was  desired  because  of  its 
simplicity. 

One  important  question  that  this  study  was  intended  to 
help  answer  is  that  of  the  feasibility  of  completing  the 
transfer  (within  acceptable  insertion  error  tolerances)  by 
using  open  loop  control,  as  opposed  to  some  form  of  closed 
loop  control  (e.g.,  explicit  guidance),  which  would  probably 
be  more  accurate,  but  also  much  more  complex;  simplicity  and 
cost  constraints  again  emphasized. 

Due  to  the  time  limit  on  this  study,  only  transfers  be- 
tween circular  orbits,  both  coplanar  acid  non-coplanar,  were 
considered.  The  emphasis  was  on  accomplishment  of  two  par- 
ticular transfers,  both  of  which  are  potentially  important 
IUS  missions.  The  first  is  from  a 160  nm  parking  orbit  at 
an  inclination  of  28.5°*  to  a geosynchronous  mission  orbit. 
The  second  involves  a transfer  from  a 160  nm  parking  orbit 
inclined  at  57°*  to  a subsynchronous  (12  hour)  mission  orbit 
with  an  inclination  of  63°. 

A computer  simulation  was  created  to  target  and  evaluate 
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the  proposed  energy  management  scheme.  The  inputs  and  outputs 
of  the  simulation  are  listed  in  Chapter  II,  and  a verbal  flow 
chart  can  be  found  in  Appendix  A.  The  following  chapters 
describe  the  formulation  of  the  guidance  scheme  which  led  to 
the  computer  simulation.  The  actual  computer  code  listing  is 
found  in  Appendix  F,  and  includes  comment  cards  highlighting 
each  important  computation. 
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II.  Formulation  of  the  Scheme 


Overview 

The  energy  management  scheme  proposed  herein  was  conceived 
with  simplicity,  practicality,  and  cost  as  the  major  considera- 
tions. In  this  scheme  the  entire  burden  of  managing  the  excess 
fuel  and  effecting  the  transfer  is  put  on  pre launch  targeting, 
so  that  during  the  bums  the  thrust  can  be  directed  along  a 
precomputed  direction  using  constant  attitude  maneuvers  only. 

The  highlight  of  this  scheme  is  its  simple  "open  loop"  de- 
sign, suggesting  minimal  on-board  equipment  for  its  execution. 
Only  six  mission  parameter  values  (outputs  of  the  targeting) 
need  be  stored  on-board  the  IUS  for  execution  of  the  transfer 
maneuver.  They  are  the  two  ignition  times,  and  two  thrust  di- 
rection angles  for  each  bum.  These  six  values  will  be  referred 
to  as  the  control  parameters!  and,  as  such,  they  could  be  im- 
plemented by  the  on-board  guidance  system  to  drive  the  IUS 
state  vector  to  match  that  of  the  mission  orbit.  Targeting, 
to  determine  the  values  of  the  control  parameters,  is  explained 
fully  in  Chapter  IV. 

Motivation  for  the  use  of  constant  attitude  thrusting 
is  due  to  the  fact  that  the  IUS  is  to  have  an  Inertial  Measure- 
ment Unit  (IMU)  which  uses  "strapdown"  gyros.  A possible  dis- 
advantage of  the  two  aforementioned  Burner  II  and  Trident  en- 
ergy management  proposals  (and  consequently  a possible  advan- 
tage of  constant  attitude  thrusting)  is  that,  in  those  schemes, 
vehicle  turning  rates  during  thrusting  can  become  quite  high, 
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and  may  adversely  affect  the  accuracy  of  the  strapdown  IMU. 
Turning  rates  of  as  much  as  7°/sec  for  the  Burner  II  (through 
a total  attitude  change  of  ±90°  during  the  burn)  were  indicat- 
ed by  the  Boeing  proposal  (Ref  2i319).  Similarily  for  the 
Trident  scheme,  if  the  vehicle  capability  is  25#  in  excess  of 
the  velocity  change  required,  then  the  vehicle  must  rotate 
through  an  attitude  change  of  125°  during  the  burn.  This  would 
produce  peak  turning  rates  of  3-^°/sec.  This  is  cited  as  a 
"significant  disadvantage  . . . which  may  affect  navigation 
accuracy  or  computation  rates  associated  with  the  strapdown 
IMU"  (Ref  3«10).  Constant  attitude  thrusting  completely  elim- 
inates this  possible  source  of  trouble. 

In  addition  to  the  decision  to  use  constant  attitude 
thrusting  only,  several  other  considerations  were  important 
in  the  early  formulation  of  this  scheme i 

1.  The  first  was  the  criterion  for  initiation  of  the 
second  stage  burnj  in  the  context  of  open  loop  control.  That 
is,  whether  the  second  bum  should  occur  at  a certain  prepro- 
grammed time,  or  at  a certain  position,  as  indicated  by  the 
on-board  navigation  equ.pment. 

As  indicated  above,  the  decision  was  made  to  base  it 
on  a predetermined  time,  so  that  the  six  control  parameters 
consist  of  the  two  start  bum  times,  and  two  thrust  angles 
for  each  bum.  The  reason  for  this  choice  is  explained  in  a 
later  section  of  this  chapter. 

2.  The  next  consideration  was  the  necessity  to  include 
finite  bum  dynamics  to  realistically  test  the  feasibility  of 
the  open  loop  design,  and  to  obtain  control  parameter  values 
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consistent  with  real  hardware. 

3*  Lastly,  since  the  IUS  must  have  the  capability  of 
accounting  for  any  mission  delays,  the  scheme  had  to  have  a 
contingency  retargeting  capability.  This  would  allow  the  trans- 
fer to  be  performed  on  consecutive  opportunities. 

The  first  idea  considered  was  to  formulate  the  problem 
using  optimal  control  theory,  where  numerical  solution  of  the 
associated  Two  Point  Boundary  Value  Problem  (TPBVP)  would  have 
automatically  accounted  for  the  finite  bums.  This  formula- 
tion would  have  selected  the  transfer  which  effectively  used 
all  the  velocity  change  capabilities  of  each  stage,  and  also 
gave  the  minimum  insertion  error  (cost  function  being  posi- 
tion and  velocity  insertion  errors).  However,  the  optimal 
control  approach  was  dropped  for  a more  flexible  procedure 
that  would  lend  more  insight  into  the  actual  maneuver  execu- 
tion, and  also  yield  more  output  information.  The  concept  of 
attempting  to  find  the  transfer  that  would  be  optimal  in  some 
sense  was  kept,  however. 

Instead  of  the  optimal  control  approach,  the  mission  con- 
straints (values  which  define  the  mission  orbit)  are  expressed 
as  nonlinear  functions  of  the  control  parameters,  and  a non- 
linear equation  solving  routine  is  used  to  search  out  the  val- 
ues of  the  control  parameters  which  cause  the  IUS  state  vector 
to  exactly  match  that  of  the  required  mission  orbit  after  ter- 
mination of  the  second  bum  (insertion). 

This  process  is  the  very  heart  of  the  scheme  developed 
in  this  study.  The  nonlinear  equation  solving  routine  accom- 
plishes exact  targeting  for  the  finite  bum  dynamics.  It  has 
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the  most  important  advantage  that  there  are  no  guidance  algo- 
rithm-generated insertion  errors  (within  the  framework  of  the 
dynamic  model). 

The  nonlinear  equation  solver  is  a general  purpose  sub- 
routine developed  by  the  mathematician  M.J.D.  Powell  (Ref  7). 

Its  callname  is  NS01A,  and  it  will  be  referred  to  here  by  that 
name . 

Early  in  the  study,  it  was  apparent  that  for  most  combi- 
nations of  velocity  change  capabilities  (AV^  and  AV2),  a range 
of  possible  transfer  trajectories  exists;  where  any  trajectory 
in  that  range  can  be  made  to  satisfy  the  energy  management  re- 
quirement. The  reason  for  this  is  that  in  most  cases  the  num- 
ber of  control  parameters  exceeds  the  number  of  mission  orbit- 
al elements  which  must  be  satisfied,  thus  introducing  extra 
degrees  of  freedom. 

The  observation  that  there  will  normally  be  a variety  of 
trajectories  available  within  a certain  range  led  back  to  the 
idea  of  selecting  an  optimal  trajectory.  The  parameter  used 
to  define  this  range  was  chosen  to  be  the  span  between  the 
minimum  and  maximum  amounts  of  plane  change  that  could  be  ac- 
complished by  the  first  bum,  and  still  satisfy  the  constraints. 

As  an  example,  for  the  geosynchronous  mission  up  to  10° 
of  the  total  plane  change  may  be  accomplished  by  the  first 
bum.  This  is  based  on  Burner  II  specifications,  which  pro- 
duce a AV^  = 9^53  ft/sec,  and  a AV2  = 7070  ft/sec  for  a 3000 
lb  payload.  So  by  sampling  this  range  at  one  degree  intervals, 
eleven  possible  transfer  trajectories  are  available  for  direct 
comparison. 
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This  is  the  standard  method  used  in  this  study  to  compare 
the  results  of  targeting  for  any  given  combination  of  energy 
capabilities,  and  AV2>  Sensitivities  to  error  inputs  are 
computed  for  each  trajectory  in  the  range,  so  that  a compari- 
son can  be  made  to  determine  which  trajectory  gives  the  smal- 
lest insertion  errors,  minimum  transfer  time,  or  whatever  the 
optimal  criterion  for  any  particular  mission  might  be. 

This  section  was  intended  to  give  the  general  reader  some 
background  and  insight  into  the  ideas  involved  in  the  formula- 
tion of  this  scheme.  All  of  these  concepts  are  explained  in 
detail  in  the  following  chapters. 

Dynamic  Model 

The  system  consisting  of  the  IUS  vehicle  undergoing  an 
orbital  transfer  about  the  earth,  is  modeled  under  the  fol- 
lowing assumptionsi 

1.  Only  two-body  equations  of  motion  apply,  with  thrust 
as  the  only  perturbative  acceleration.  That  is,  any  perturba- 
tions due  to  solar  radiation  pressure,  and  the  gravitational 
effects  of  the  sun,  moon,  and  other  celestial  bodies  are  as- 
sumed negligible.  The  restricted  two-body  equations  of  motion 
are  presented  in  Appendix  C. 

2.  An  inverse  square  gravitational  field  applies  about 
a spherical  earth  (i.e.,  earth  oblateness  effects  are  negligi- 
ble) . 

3.  The  mass  flow  rate  (bum  rate)  of  each  solid  rocket 
motor  is  assumed  constant  with  time,  thus  producing  a constant 
thrust.  Initial  thrust  buildup  and  final  thrust  tail-off  ef- 
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fects  are  assumed  insignificant. 

4.  When  computing  the  performance  characteristics  of 
each  stage  needed  for  the  impulsive  targeting  first  approxi- 
mation, the  ideal  velocity  equation  (with  later  corrections 
for  finite  bum  losses)  is  assumed  to  apply « 


AV  = I 


ft) 


(2-1) 


where  AV  is  the  ideal  velocity  capability  of  the  stage ; I 

sp 

is  the  specific  impulse;  gQ  is  the  gravitational  constant; 
mQ  is  the  mass  prior  to  ignition;  and  mf  is  the  final  mass 
after  burnout. 

Although  the  computer  simulation  used  for  this  study  (and 
consequently  the  method  itself)  will  accept  any  IUS  vehicle 
specifications,  the  Burner  II  values,  as  given  in  Reference  2, 
are  used  throughout  to  standardize  the  results. 

The  assumptions  and  constraints,  under  which  an  error 
analysis  of  this  scheme  is  accomplished,  are  described  in 
the  next  section. 


Error 


An  underlying  objective  of  this  study  is  to  determine  the 
feasibility  of  completing  the  orbital  transfer  using  simple 
open  loop  control;  under  the  presumption  that  any  additional 
software  or  Reaction  Control  System  (RCS)  correction  bums, 
may  be  unnecessary.  If  insertion  errors  could  be  kept  within 
an  acceptable  range  (by  using  the  most  optimal  trajectory), 
then  implementation  of  this  type  of  a scheme  might  prove  fea- 
sible . 
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A major  point  of  emphasis  concerning  this  particular 
scheme  is  that  within  the  framework  of  the  model  just  stated, 
it  is  exact.  That  is,  if  there  were  none  of  the  below  listed 
unmodeled  disturbance  inputs,  there  would  be  no  insertion  er- 
ror after  execution  of  the  scheme. 

External  disturbances  will  be  present  to  some  degree,  how- 
ever, and  will  introduce  errors  into  the  transfer  maneuver. 

The  main  error  sources  are  as  follows: 

1.  IMU  errors  (Ref  6) 

a.  Initial  alignment  errors 

b.  Gyro  drift-rate  bias 

c.  Acceleration-sensitive  gyro  drift 

d.  Accelerometer  bias 

e.  Accelerometer  scale  factor 

f.  Gyro  torquer  scale  factor 

g.  Gyro  input  axis  alignment 

h.  Gyro  torquer  asymmetry 

2.  Velocity  change  perturbations 

a.  Vehicle  structure  and  fuel  weight  deviations 

b.  Specific  impulse  (I__)  deviations 

sp 

c.  Thrust  profile  fluxuations 

3.  Gravity  perturbations 

a.  n-Body  disturbances 

b.  Earth  oblateness  effects 

4.  Solar  radiation  pressure 

Due  to  the  time  constraint  involved  in  this  study,  all  of 
the  above  error  sources  could  not  be  included  in  the  accuracy 
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analysis  of  the  scheme.  Those  selected  can  be  neatly  summa- 
rized as  thrust  vector  errors*  and  constitute  the  most  signif- 
icant disturbances.  Thrust  vector  errors  originate  from  sourc- 
es 1.  and  2.c. 

The  overall  effect  of  IMU  errors  is  to  cause  an  error  in 
the  direction  of  the  applied  thrust.  This  may  be  termed  thrust 
misalignment  error.  Error  source  2.c.  arises  from  variations 
in  the  mass  flow  rate  (bum  rate)  of  the  engines.  This  causes 
an  error  in  the  thrust  magnitude,  which  in  turn  perturbs  the 
velocity  change  acceleration  profile,  ultimately  causing  an 
error  in  the  burnout  position. 

In  summary,  the  two  general  sources  cf  error  considered 
to  affect  this  scheme  are  deviations  in  the  thrust  vector  mag- 
nitude and  direction.  Insertion  error  sensitivities  due  to 
both  thrust  misalignment  and  thrust  magnitude  deviations  are 
computed  for  each  trajectory  targeted. 

Choice  of  a Transfer  Scenario 

Once  it  was  decided  that  guidance  would  be  performed  us- 
ing constant  attitude  bums,  it  was  necessary  to  decide  how 
best  this  could  be  implemented  in  hardware.  The  choice  of 
selecting  the  control  parameters,  based  on  two  bum  times, 
appeared  to  be  the  best  method,  as  explained  here. 

Assuming  that  the  first  stage  ignition  will  occur  at  the 
proper  position  in  any  scenario  (due  to  the  proximity  of  the 
IUS  to  the  orbiter  vehicle  with  its  position  well  known),  then 
there  are  two  different  possibilities  for  initiation  of  the 
second  bum.  One  where  second  stage  ignition  occurs  when  the 
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on-board  IMU  indicates  it  should,  or  one  where  second  stage 
ignition  occurs  according  to  an  on-board  clock  at  a certain 
preprogrammed  time.  The  question  obviously  was  which  scenar- 
io would  be  the  most  accurate.  That  is,  which  case  would  be 
least  sensitive  to  thrust  vector  errors. 

To  investigate  this  question  two  similar  computer  sim- 
ulations were  developed.  In  both  cases  Hohmann  transfer  ve- 
locity change  capabilities  were  assumed  and  the  impulsive  ap- 
proximation was  used.  An  ideal  IMU  was  assumed  (with  no  drift, 
etc.),  and  an  alignment  error  present  in  both  cases. 

The  first  simulation  initiated  the  second  burn  when  ei- 
ther the  IMU  indicated  that  the  proper  altitude  had  been 
reached,  or  that  180°  of  transfer  angle  had  been  completed. 

The  second  simulation  precomputed  the  transfer  time  of  flight 
and  initiated  the  second  burn  at  that  instant  along  the  trans- 
fer trajectory.  Insertion  error  sensitivities  were  computed 
for  each  case  and  combined  into  values  for  position  insertion 
error  and  velocity  insertion  error.  The  transfers  tested  were 
between  a 160  nm  parking  orbit  and  a synchronous  orbit  (at 
19*323  nm) . Plane  changes  of  0°,  28. 5°  and  57°  were  accom- 
plished. In  all  cases  the  TOF  initiated  second  bum  performed 
more  accurately,  by  about  a factor  of  two,  as  shown  in  Tables 
I and  II  for  a one  milliradian  misalignment  of  the  IMU  axes 
during  both  bums. 

The  logical  choice  of  a transfer  scenario  then,  from  these 
results,  was  one  based  on  transfer  time  of  flight. 

gxp.lanaU.9n  al  ifcg.  fae launch  Targeting 
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Table  I. 


Insertion  Errors  with  IMU 
Initiated  Second  Bum 


Plane 

Change 

Position 
Error  (nm) 

Velocity 
Error  (ft/sec) 

0° 

22.2 

15.7 

28.5° 

22.2 

14.5 

57° 

22.2 

11.1 

Table  II. 

Insertion  Errors  with  TOF 
Initiated  Second  Bum 


Plane 

Change 

Position 
Error  (nm) 

Velocity 
Error  (ft/sec) 

0° 

12.6 

7.5 

28.5° 

12.6 

6.3 

57° 

12.6 

2.9 

The  final  form  for  this  scheme  followed  directly  from  the 
choice  of  a transfer  scenario  based  on  time  of  flight  as  the 
criteria  for  initiating  the  second  bum.  Certain  conditions 
needed  to  be  defined,  however,  in  order  to  lay  the  framework 
in  which  that  scenario  could  be  executed.  These  conditions 
become  the  inputs  to  the  computer  simulation  used  to  target 
the  transfer.  Both  the  inputs  and  outputs  of  the  targeting 
program  were  chosen  to  be  as  follows i 

Inputs 

1.  IUS  vehicle  specifications  (given  by  stage) 
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is,  the  values  of  t^  and  t^2  are  given  in  the  number 
of  seconds  past  epoch. 

3.  The  thrust  direction  angles  are  referenced  to  the 
geocentric-equatorial  inertial  frame. 

4.  Only  simple  plane  changes  are  accomplished.  This  re- 
quires that  either  and  Q2  are  equal,  or  that  one 
or  both  are  undefined  (equatorial  orbit),  so  that 
all  the  required  plane  change  is  just  equal  to  the 
difference  in  the  inclinations. 

5.  If  the  parking  orbit  and  mission  orbit  are  non-co- 
planar.  then  the  targeting  is  accomplished  to  place 
the  IUS  in  the  specified  mission  orbit  only  (i.e., 
the  point  of  insertion  is  not  constrained).  If  the 
orbits  are  coolanar.  then  the  targeting  automatical- 
ly accomplishes  a rendezvous  between  the  IUS  and 
the  target  position  in  the  mission  orbit. 

Properly  defining  the  epoch  time  (tQ)  allows  the  orbital 
data  to  be  expressed  in  the  form  of  the  inputs  above.  The 
scheme  presupposes  that  the  parking  orbit  is  already  estab- 
lished, and  the  position  of  the  IUS  in  that  orbit  is  accurate- 
ly known.  Thereafter,  the  epoch  time  is  defined  to  be  any  one 
of  the  times  (the  particular  one  chosen  by  the  user)  when  the 
IUS  crosses  the  line  of  the  ascending  node  while  in  the  park- 
ing orbit.  If  the  parking  orbit  is  equatorial,  then  epoch 
becomes  the  time  that  the  IUS  is  positioned  along  the  X axis 
of  the  geocentric-equatorial  frame.  With  this  definition  of 
epoch,  the  true  longitude  of  the  IUS,  in  the  parking  orbit 
(f01),  is  always  equal  to  Thus,  only  need  be  specified 
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to  fully  fix  the  IUS  position  as  a function  of  time  in  that 
orbit. 

The  perifocal  coordinate  system  (PQW  frame)  of  either  the 
parking  orbit  or  mission  orbit  as  used  here  will  have  the  P 
axis  pointing  along  the  line  of  the  ascending  node,  if  it 
exists,  or  if  the  orbit  is  equatorial,  along  the  X axis.  These 
frames  and  angular  relationships  are  illustrated  in  Figure  2-1 


as  they  might  be  at  some  epoch  time  (tQ),  for  a geosynchronous 
mission.  The  vector  ric(t)  tracks  the  IUS  in  the  parking  or- 
bit, and  the  vector  r2C(t)  tracks  the  target  position  in  the 
mission  orbit  (when  rendezvous  is  to  be  accomplished). 


It  is  an  important  characteristic  of  this  scheme  that  a 
value  for  *02  (i.e.,  an  epoch  time,  and  target  position)  al- 
ways be  specified.  This  is  necessary  so  that  the  six  control 
parameters  will  have  a reference  to  which  they  can  be  related. 
If  the  insertion  point  of  the  IUS  into  the  mission  orbit  is 
unconstrained  (for  a coplanar  transfer),  or  if  the  transfer 
is  non- coplanar,  then  an  arbitrary  value  for  *02  may  be  used 
(for  instance,  iQ2  = 0). 

It  happens  that  transfer  opportunities  are  regularly  re- 
petitive between  non-coplanar  orbits  (when  rendezvous  is  un- 
necessary), and  synodically  repetitive  between  coplanar  or- 
bits (to  satisfy  a rendezvous).  This  allows  contingency  re- 
targeting to  be  accomplished,  since  an  endless  list  of  pos- 
sible mission  start  times  are  available.  Thus,  the  scheme  is 
programmed  to  generate  the  values  of  the  control  parameters 
associated  with  the  first  five  mission  opportunity  times  af- 
ter epoch,  for  each  transfer  trajectory  targeted.  The  reason 
for  doing  this  is  that,  if  a mission  delay  were  to  occur,  the 
control  parameters  could  be  automatically  reset  to  the  next 
sequential  set  of  values. 

The  algorithm  used  to  accomplish  this  targeting  task  is 
presented  in  the  next  two  chapters. 
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III.  Energy  Management  by  Trajectory  Matching 


Various  possibilities  have  seen  suggested  for  managing 
the  excess  energy  of  the  fixed-impulse  solid  rocket  motors. 

The  prominent  ideas  have  included  propellant  offloading,  ad- 
ding ballast,  and  attitude  modulation  of  the  thrust  vector* 
all  of  which  have  major  drawbacks.  This  chapter  explains  the 
energy  management  technique  employed  by  this  scheme,  which  in- 
volves selecting  a non-Hohmann  transfer  trajectory  which  "match 
es"  the  fixed  energy  capabilities  (AV^  and  AV2)  of  the  IUS  ve- 
hicle. 

For  a AV^^  and  AV2  combination,  both  in  excess  of  the 
Hohmann  (minimum  energy)  values  for  any  particular  mission, 
a non-Hohmann  transfer  is  usually  possible.  A simple  copla- 
nar  non-Hohmann  transfer  is  illustrated  in  Figure  3-1* 

For  any  orbital  transfer  maneuver,  the  insertion  con- 
ditions may  be  related  to  the  conditions  at  the  first  burn  by 
using  the  impulsive  velocity  change  approximation,  and  a ser- 
ies of  orbital  transfer  equations.  This  series  of  equations 
is  contained  in  Appendix  B as  Equations  (B-ll)  through  (B-28) . 
Within  these  equations  a relationship  exists  between  possible 
values  of  AV^  and  AV2  in  the  form  of 

AV2(AVlt  virlt  ij;2)  (3-1) 

where  for  any  value  of  AV^,  ^ is  the  associated  first  burn 
flight  path  angle,  and  t|J2  is  the  plane  change  angle  accom- 
plished by  the  first  burn.  The  consequent  AV2  is  the  value 
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Figure  3-1.  Coplanar  non-Hohmann  Transfer 


of  the  second  bum  velocity  increment  which  completes  any  re- 
maining plane  change  and  causes  insertion  into  the  mission  or- 
bit. 

In  order  to  determine  which  combinations  of  AV.^  and  AV2 
would  allow  trajectory  matching  as  a means  of  energy  manage- 


ment, plots  were  made  of 


AV  combinations  for  each  of 


the  reference  missions  listed  in  Chapter  I.  These  plots  are 
shown  in  Figures  3-2,  3-3»  and  3-4,  for  the  geosynchronous, 
subsynchronous,  and  geosynchronous  coplanar  transfers,  respec- 


tively. 


Each  of  these  plots  was  obtained  by  accomplishing  a two- 
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dimensional  mapping  through  all  allowable  combinations  of  tjj^ 
and  for  various  values  of  AV^  using  the  orbital  transfer 
relationships  of  Equations  (B-ll)  through  (B-28) . Using  this 
procedure,  the  boundaries  of  the  allowable  regions  become  the 
maximum  and  minimum  values  that  AV2  may  take  on,  for  a set 
value  of  AV^ . 

It  happens  that  the  maximum  possible  value  for  AV2  (giv- 
en the  value  of  AV^)  occurs  where  both  ijj-^  and  iJj2  are  equal  to 
zero.  The  minimum  possible  AV2  occurs  where  ij;2  has  its  max- 
imum value,  and  ^ (second  burn  flight  path  angle)  becomes 
equal  to  zero.  A trajectory  is  possible  for  any  value  of  AV2 
between  these  extremes. 

Superimposed  on  Figures  3-2  and  3-3  are  plots  of  the  A)^ 
and  AV2  resulting  from  various  payload  weights,  using  the 
Burner  II  as  the  IUS  vehicle,  with  a first  stage  fuel  load  of 
17,300  lbs  and  a second  stage  fuel  load  of  4?00  lbs.  It  is 
apparent  that  the  excess  energy  balance  between  AV1  and  AV2 
is  primarily  dependent  on  the  payload  weight. 

The  significance  of  these  plots  is  that  transfer  trajec- 
tories may  always  be  found  that  match  any  combination  of  AV's 
selected  from  within  the  allowable  regions.  In  the  case  of 
a coplanar  transfer,  the  functional  relationship  of  (3-1)  re- 
duces to 

AV2(^1)  (3~2) 

for  a fixed  value  of  AV^.  The  trajectory  matching  procedure 
then  consists  of  finding  the  solution  of 
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(3-3) 


= av2* 

£ 

where  AV2  is  the  fixed  (known)  value. 

The  analytical  relationship  of  contained  in  Equa- 

tions (B-ll)  through  (B-28),  is  highly  nonlinear  and  transcen- 
dental in  nature,  and  must  be  solved  iteratively.  To  illus- 
trate this,  a closed  form  expression  for  AV2(\Jj^)  is  derived  in 
Appendix  B,  for  a coplanar  transfer. 

For  a transfer  between  coplanar  orbits,  only  one  trajec- 
tory exists  which  matches  a given  AV1  and  AV2  combination,  and 
remains  within  the  plane  of  the  orbits.  It  is  often  possible, 
however,  to  gain  flexibility  by  thrusting  out  of  plane  with  a 
portion  of  AV1  and  back  into  plane  with  a portion  of  AV2>  It 
is  interesting  to  note,  though,  that  the  minimum  time  of  flight 
trajectory  is  always  the  one  which  remains  in  plane.  The  rea- 
son for  this  is  that,  in  this  case,  all  the  available  energy 
remains  in  the  plane  of  the  transfer  trajectory. 

For  a transfer  between  non- coplanar  orbits,  an  additional 
degree  of  freedom  is  introduced  through  the  addition  of  4r2,  as 
expressed  in  Equation  (3-1)*  A particular  combination  of  AV^ 
and  AVg  may  now  yield  a variety  of  different  trajectories,  de- 
pending on  the  amounts  of  plane  change  accomplished  during  each 
burn.  In  this  case,  there  will  be  a range  of  values  for  ilr2, 
any  one  of  which  will  allow  a solution  to  the  remaining  rela- 
tionship of  Equation  (3-3)* 

The  actual  range  over  which  2 may  take  on  values  depends 
on  the  excess  energy  balance  between  AV^  and  AV2  (i.e.,  their 
coordinates  within  the  allowable  region) . The  upper  and  lower 
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limits,  on  the  range  of  values  i|/2  may  have,  are  always  given 
byi 


0° 


S ^2  * ^2max 


(3-4) 


A relationship  for  calculating  is  given  in  Appendix 

B.  The  actual  range  of  iJj2  (for  finite  burn  trajectories,  etc.) 
will  always  be  somewhat  less  than  this. 

Thus,  trajectory  matching  using  the  impulsive  velocity 
change  approximation  becomes  an  important  early  step  in  the 
targeting  process.  Acceptable  trajectories  are  found  by  fix- 
ing a value  of  i(r2,  which  is  the  parameter  used  to  define  the 
range  of  usable  trajectories,  and  then  accomplishing  a one~ 
dimensional  search  on  AV2,  for  varying  values  of  ij^.  When 
the  solution  to  Equation  (3-3)  has  been  bracketed  by  this 
search  procedure,  a Regula-Falsi  (Ref  4«178)  iteration  is  ac- 
complished to  refine  it  exactly.  During  the  targeting  compu- 
tations (described  in  Chapter  IV)  each  value  of  iJ/2,  within  the 
range  (usually  at  1°  intervals),  is  sampled  to  generate  all 
the  possible  trajectories. 

The  amount  of  flexibility  in  choosing  possible  trajec- 
tories depends  entirely  on  the  range  of  values  tjj2  may  take 
on,  which  in  turn  depends  on  the  energy  balance  between  AV^ 
and  AV2.  Only  one  trajectory  is  possible  for  a AV  combina- 
tion lying  on  a boundary  of  the  allowable  regions.  This  be- 
comes the  case  for  the  maximum  payload  combination. 

r 
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IV.  Targeting 


Overall  Objective 

The  mission  objective  is  to  complete  an  orbital  trans- 
fer between  a low  altitude  pai'king  orbit  and  a higher  alti- 
tude mission  orbit.  The  mission  orbit  can  be  defined  by 
certain  constraints#  or  "hit  conditions",  to  which  the  IUS 
state  vector  must  be  driven  in  order  to  accomplish  inser- 
tion into  that  orbit.  The  form  of  the  mission  constraints 
used  in  this  targeting  algorithm  are  chosen  to  conform  with 
the  energy  management  technique  described  in  the  last  chap- 
ter. 

For  a vehicle  with  fixed-impulse  capability,  the  pay- 
load  weight  is  the  governing  parameter  in  fixing  the  actual 
values  of  AV1  and  AV2,  and  the  consequent  excess  energy  bal- 
ance. With  excess  energy,  a variety  of  trajectories  are  pos 
sible  which  would  satisfy  the  mission  constraints.  The  pa- 
rameter \Jj2  (plane  change  angle  accomplished  by  the  first 
bum)  was  chosen  to  define  the  range  of  transfers  possible 
using  the  available  AV1  and  AVg.  The  idea  was  to  fix  values 
of  ijr2  at  one  degree  intervals  through  the  range  so  that  a 
comparison  of  the  resulting  transfer  trajectories  could  be 
made. 

Specifying  \\f2  for  a given  transfer  is  equivalent  to  con 
straining  the  inclination  of  the  transfer  trajectory  (iT). 
For  a fixed  amount  of  first  bum  plane  change  (i^2)  • 
inclination  of  the  parking  orbit  (i^)  specified,  the  incli- 
nation of  the  transfer  orbit  is  the  fixed  sum  of  these  two 
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angles. 

In  the  targeting  algorithm,  the  IUS  vehicle  state  vec- 
tor (three  components  of  position,  and  three  components  of 
velocity)  becomes  a function  of  the  six  control  parameters. 
Where  necessary,  the  components  of  the  IUS  state  vector  are 
converted  into  the  equivalent  orbital  elements  for  direct  com- 
parison with  the  mission  constraints.  The  function  of  the 
targeting  algorithm  is  to  generate  the  values  of  the  control 
parameters  which  drive  the  IUS  state  vector  (and  its  cor- 
responding orbital  elements)  to  match  the  values  of  the  mis- 
sion constraints.  The  algorithm,  as  coded,  will  only  tar- 
get transfers  between  circular  orbits.  It  would  be  straight- 
forward, however,  to  include  coding  which  would  handle  el- 
liptical orbits  also. 

Mission  Constraints  and  Variables 

Six  constraints  are  required  to  completely  specify  an 
orbit,  in  the  most  general  case.  Five  of  the  constraints 
may  be  constants  which  specify  the  size,  shape,  and  orien- 
tation of  the  orbit.  The  sixth  constraint  specifies  the 
position  within  that  orbit,  and  is,  therefore,  a function 
of  time.  The  classical  orbital  elements  are  the  parameters 
most  often  used  to  specify  an  orbit.  The  mission  constraints 
in  this  algorithm,  are  a modification  of  the  classical  orbit- 
al elements.  The  IUS  vehicle  must  satisfy  some  or  all  of 
these  constraints  after  termination  of  the  second  stage 
bum  (depending  on  the  type  of  transfer  - as  will  be  ex- 
plained) . 
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Mission  Constraintsi  (hit  conditions) 


f 


1. 

r2 

- magnitude  of  the  mission  orbit  position 
vector 

2. 

V2 

- magnitude  of  the  mission  orbit  velocity 
vector 

3- 

e2 

- eccentricity  of  the  mission  orbit 

4. 

i2 

- inclination  of  the  mission  orbit 

5- 

°2 

- longitude  of  the  ascending  node  of  the  mis- 
sion orbit 

6. 

» 

TA 

- transfer  angle 

Of  the  six  mission  constraints  shown*  the  first  five 
are  constants  (both  r2  and  V2  are  constant  for  a circular 
orbit),  but  the  sixth  is  a function  of  time  and  needs  some 
explanation.  In  essence,  TA(t)  is  the  selected  form  of  the 
orbital  element  which  insures  the  insertion  into  the  mis- 
sion orbit  at  the  desired  position  in  that  orbit.  As  used 
here,  TA(t)  is  defined  to  be  the  angle  between  the  two  vec- 
tors  r1(.(t)  and  r^ft  + TOF  ),  where  "t"  is  any  particular 
time,  and  TOF  is  the  orbital  transfer  timej  that  is,  the 
time  from  initiation  of  the  first  stage  bum  to  termina- 
tion of  the  second  stage  bum. 

Before  the  transfer,  the  position  of  the  IUS  vehicle 
is  tracked  by  the  vector  ric(t),  and  the  target  position 
in  the  mission  orbit  is  tracked  by  the  vector  r2C(t).  Dur- 
ing the  actual  transfer,  the  target  position  will  have  shift- 

_ _ # 

ed  along  its  orbit  from  r2C(ttol)  to  **20  (t^  + TOF  ),  where 
t^  is  the  time  the  transfer  began  (i.e.,  first  stage  igni- 
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tion  time).  This  is  illustrated  in  Figure  4-1  for  the  case 
of  a simple  coplanar  transfer. 

Satisfaction  of  the  sixth  constraint  really  means  find- 
ing the  correct  first  stage  ignition  time,  such  that 

T A ( t^i ) = TA*  (4-1) 

# * — _ 
where  TA  is  the  angle  between  the  vectors  and  r2C 

(tbi  + TOF*).  Thus,  it  is  apparent  that  the  constraint  TA* 

is  dependent  for  its  actual  value  on  the  transfer  time  of 

# 

flight  TOF  , which  in  turn  is  dependent  on  the  particular 

transfer  trajectory. 

* 

Since  TOF  is  really  dependent  on  the  two  burn  times, 

# # 

the  true  functional  form  of  TA  really  is  TA  (t^,  t^)  • 


TOF* 


Figure  4-1.  Transfer  Angle 
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Satisfaction  of  TA  may  be  thought  of  as  insuring  the  proper 
"phase  angle’*  at  departure!  to  effect  a rendezvous  with  the 
target  position  in  the  mission  orbit. 

The  six  control  parameters  are  the  variables  used  to 
satisfy  the  mission  constraints#  and  are  listed  herei 

Control  Parameters! 

1.  ttl  - time  of  first  stage  ignition  (referenced  to 
epoch,  t0  = 0) 

rpi  - l inertial  thrust  direction  angles  of  the 

2 ^ | first  stage  bum 

4.  tt2  - second  stage  ignition  time  (a  function  of 

transfer  TOF) 

5*  $3  ) inertial  thrust  direction  angles  of  the 

6 - jsecond  stage  bum 

4 

The  thrust  direction  angles  specify  the  constant  atti- 
tude inertial  orientations  of  the  thrust  vectors  of  each 
stage,  T^  and  T2.  The  angles  cp1  and  cp2  define  the  direc- 
tion of  T1§  as  shown  in  Figure  4-2.  Similarly,  the  angles 
qp^  and  $4  specify  the  orientation  of  Tg. 
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The  important  point  in  the  relationship  between  the  mis- 
sion constraints  and  the  control  variables  is  that  unique 
values  for  the  six  control  parameters  specify  unique  values 
for  the  six  constraints.  That  is.  any  combination  of  values, 
even  selected  arbitrarily,  for  the  control  parameters  will 
produce  some  unique  orbit  (not  necessarily  circular)  after 
the  second  bum.  The  task  of  the  targeting  algorithm  then 
becomes  one  of  generating  the  proper  values  of  the  control 
parameters  such  that  their  values  drive  the  IUS  vehicle  to 
produce  an  orbit  after  second  stage  burnout  which  satisfies 
the  constraints. 

There  are  many  ways  to  attack  this  problem,  but  the 
method  decided  upon  here  (which  attempts  to  be  conceptually 
straightforward),  is  to  use  two  distinct  steps  in  the  solu- 
tion process.  The  first  step  is  to  simply  target  the  trans- 
fer in  a general  form,  where  only  four  of  the  constraints 
are  satisfied.  The  four  to  be  satisfied  in  this  first  step 
of  targeting  are  r2,  V2,  e2,  i2»  with  n2  and  TA  left  free. 

The  second  step  essentially  satisfies  proper  phasing. 

The  decision  to  target  the  transfer  in  two  steps  can  be 
explained  best  by  looking  at  the  functional  relationship  be- 
tween the  constraints  and  variables! 

1.  r2^l*  $2*  ^b2'  ^3*  ^P4^ 

2.  ^2^1*  ^2*  ^b2*  ^3*  *£4) 

3*  e2^1*  $2'  ^b2*  ^3*  ^4) 

4.  i 2 (cp]_ * cpg#  ^)2*  Vy  ^4^ 

5*  ^2^^bl'  ^1’  ^2*  ^b2*  ^3*  ^4^ 

6.  TA(t^^,  t^jg) 
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Thus,  it  is  apparent  that  there  is  decoupling  between 
the  first  four  constraints  and  the  last  two  constraints, 
with  respect  to  the  variable  t^.  The  first  four  constraints 
are  not  a function  of  t^.  In  other  words,  they  are  inde- 
pendent of  the  position  in  the  parking  orbit  at  which  first 
stage  ignition  occurs  and  the  transfer  is  initiated.  Ver- 
ification of  the  above  constraint  variable  relationships 
through  the  interrelating  equations  is  straightforward  and 
will  not  be  shown  here. 

Since  there  are  six  constraints  and  an  equal  number  of 
variables,  it  might  to  possible  to  find  a unique  solution 
for  the  variables  which  would  satisfy  all  six  constraints, 
for  even  the  most  general  non-coplanar  transfers.  However, 
due  to  the  nature  of  the  initial  conditions  (Chapter  II, 
last  section)  and  the  energy  management  requirements  incor- 
porated in  the  governing  relationships,  a straightforward 
way  of  solving  for  values  of  the  control  parameters  which 
would  satisfy  all  six  constraints  in  every  case,  was  not 
apparent.  From  the  familiarization  gained  in  this  study  it 
seems  most  likely  that  a unique  solution  is  only  possible 
under  certain  circumstances. 

The  reason  for  the  difficulty  in  finding  a unique  so- 
lution to  satisfy  all  six  constraints,  is  that  one  degree 
of  freedom  is  "lost"  when  the  parking  orbit  is  specified 
before  targeting  (as  is  always  the  case  in  this  scenario). 
Although  the  parameter  t^  is  free,  it  is  really  the  vector 
?lC(t)  which  must  be  completely  free  in  order  for  a solution 
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to  be  possible  which  satisfies  all  six  constraints.  When 
the  IUS  position  in  the  parking  orbit  is  specified  before 
targeting,  r^(t)  is  no  longer  free  with  respect  to  time. 

There  is,  however,  a certain  limited  degree  of  freedom 
generated  by  the  range  of  values  i{f2  may  have,  which  in  turn 
gives  some  flexibility  in  choosing  a trajectory  with  a cer- 
tain transfer  angle.  Here,  even  at  best,  a unique  solution 
would  involve  a very  narrow  launch  window  for  the  space  shut- 
tle orbiter  itself,  and  even  then  may  allow  only  one  mission 
opportunity,  thus  precluding  a retargeting  capability. 

Since  most  missions  will  probably  not  involve  a ren- 
dezvous requirement  anyway,  the  targeting  algorithm  devel- 
oped here  assumes  an  arbitrary  parking  orbit  already  estab- 
lished. When  satisfaction  of  the  sixth  constraint  (ren- 
dezvous) is  required  between  non-coplanar  orbits,  it  may 
still  be  accomplished  by  a station-keening  maneuver  within 
the  parking  orbit  to  obtain  proper  phasing  for  the  targeted 
trajectory. 

Thus,  the  path  chosen  in  the  formulation  of  this  scheme 
was  to  satisfy  only  five  of  the  six  constraints  (through  the 
targeting  algorithm)  for  general  non-coplanar  transfer  mis- 
sions. In  step  one  of  targeting  the  first  four  constraints 

are  satisfied,  and  in  targeting  step  two  (phasing)  either  Cl2 

# 

or  TA  is  satisfied  depending  on  the  mission.  If  the  trans- 
fer is  between  coolanar  orbits,  n2  is  never  a constraint 
since  it  is  always  satisfied,  and  TA  is  enforced  in  this 
case  as  the  fifth  constraint.  If  the  mission  is  a non-co- 
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planar  transfer,  then  02  is  satisfied  in  step  two  and  TA 
(insertion  position)  is  left  free. 

Targeting  Step  One 

The  first  step  of  the  targeting  process  involves  find- 
ing the  geometry  of  a trajectory  which  will  satisfy  the  first 
four  constraints,  without  regard  to  proper  phasing.  Leaving 
the  last  two  constraints,  02  and  TA  , free  in  this  step  leaves 
four  constraints  as  functions  of  five  variables* 

1.  r2^1'  ^2*  ^b2*  ^3’  ^4^ 

2.  ^2^1’  ^2’  ”H)2’  ^3*  ^4^ 

3*  e2^1*  ^2*  ^b2’  ^3*  ^4^ 

4.  i 2 (''Pi » ^2*  ^b2*  ^3*  ^4^ 

In  order  to  solve  for  unique  values  of  the  five  vari- 
ables which  will  satisfy  these  four  constraints,  an  addi- 
tional constraint  equation  (which  is  consistent  with  these 
four)  is  necessary.  A logical  and  very  convenient  choice 
for  the  additional  constraint  was  to  use  the  inclination 
of  the  transfer  trajectory  iT(cplt  Cp2)»  since  it  became  an 
additional  constraint  anyway  when  a value  for  4*2  was  fixed. 

In  this  sense,  iT  might  more  aptly  be  called  a design  pa- 
rameter. since  its  value  is  selected  freely  (from  within 
the  limits  of  its  range).  Once  set,  however,  iT  functions 
as  another  constraint  that  the  targeting  must  satisfy. 

Since  the  first  step  of  targeting  involves  constraints 
which  can  be  satisfied  regardless  of  where  t occurs  in 
the  parking  orbit,  a "local"  inertial  frame  (XYZ^)  is  ini- 
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tially  used.  This  local  frame  (as  it  will  be  called)  is 
referenced  to  the  first  burn  point,  as  shown  in  Figure  4-3. 


r 


( 


The  axis  always  passes  through  the  point  of  the  first 
stage  ignition.  The  frame  is  oriented  such  that  the  in- 
clination of  the  parking  orbit  is  always  zero,  and  the  in- 
clination of  the  mission  orbit  is  just  equal  to  the  abso- 
lute difference  between  the  actual  parking  and  mission 
orbit  inclinations. 

In  this  initial  step,  the  impulsive  velocity  change 
approximation  is  used  so  that  a trajectory  may  be  targeted 
(in  the  local  frame)  to  satisfy  iT,  r2,  V2,  e2,  and  i2, 
using  the  orbital  transfer  relationships  of  Equations  (B-ll) 
through  (B-28).  In  this  approximation  the  thrust  direction 


and  associated  AV  direction  are  by  definition  the  same. 
Through  these  relationships,  the  impulsive  directions  of 
AVf  (T 1 direction)  and  AV2  (T  2 direction) , as  expressed  in 
the  local  frame,  are  found  which  match  an  energy  management 
transfer  trajectory  for  a set  value  of  tj;2.  From  this,  the 
four  thrust  direction  angles  are  determined  with  respect  to 
the  local  frame . The  consequent  transfer  time  of  flight 
(TOF)  becomes  an  approximation  to  the  second  stage  bum  time 
(tb2)»  as  measured  from  a first  bum  referenced  to  tbi  = 0. 

After  approximate  values  for  cp^,  rp2,  t^2,  cp^,  and  cp^ 
are  obtained  through  impulsive  targeting,  their  actual  val- 
ues for  the  real  finite  burn  dynamics  are  found  using  the 
nonlinear  equation  solving  routine  (NS01A),  coupled  with 
numerical  integration  of  the  equations  of  motion.  Once 
NS01A  has  generated  values  for  the  five  control  parameters 
which  correspond  to  the  minimum  value  for  t|;2  (amount  of 
first  bum  plane  change),  targeting  for  subsequent  values 
of  ijr2  is  expedited  by  using  the  values  of  the  control  pa- 
rameters generated  for  the  previous  t|j2,  as  the  initial  in- 
puts to  NS01A.  This  gives  shorter  iteration  times,  as  op- 
posed to  recompleting  impulsive  targeting  at  each  step  of  2 . 

NS01A 

Since  the  five  constraint  variable  relationships  used 
in  step  one  of  targeting  are  highly  nonlinear,  a numerical 
solution  technique  must  be  employed.  NS01A,  a Fortran  sub- 
routine developed  by  Powell,  is  a highly  effective  numeri- 
cal algorithm  which  solves  a set  of  nonlinear  algebraic  equa- 
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tions.  The  NS01A  subroutine  itself  is  a standard  listing 
which  calls  on  another  subroutine  (call  named  CALFUN)  in 
which  the  particular  equations  to  be  solved  are  contained. 
Only  the  fundamental  characteristics  of  NS01A  are  described 
here.  A detailed  explanation  of  the  algorithm  and  the  ac- 
tual Fortran  listing  can  be  found  in  Reference  7. 

The  nonlinear  equations  to  be  solved  by  NS01A  are  con- 
tained in  CALFUN  (which  the  user  must  write).  The  equations 
must  be  expressed  in  the  standard  form 

F(Y)  = 0 (4-2) 


where  Y is  an  n-vector  of  the  n unknowns,  and  F denotes  the 
set  of  n nonlinear  functions,  each  expressed  as  the  differ- 
ence between  the  current  value  and  the  desired  value  of  that 
function.  Expressed  in  the  CALFUN  format  of  (4-2),  the  non- 
linear equations  to  be  solved  by  NS01A  become* 


F,  = 


Fo  = 


F,  = 


F,.  = 


F-  = 


= iTm 

- iT 

* r2(?) 

‘ r2C 

■ v7) 

‘ V2C 

= e2(Y) 

- 0 

’ i2(Y) 

" i2 

(4-3) 


The  values  on  the  right  of  the  minus  sign  are  the  actual  val- 
ues of  the  constraints  for  that  particular  transfer.  They 
are  usually  referred  to  as  the  hit  conditions.  The  Y vec- 


of  unknowns  is 


thrust  direction  angles  ex- 
pressed in  the  local  inertial 
frame 


(4-4) 


Y1  " *i 

Y2  = ^2 

Y3  = ^3 
= <P4 

= TOF  (equivalent  to  ^b2^ 

Given  the  equations  from  CALFUN  in  the  form  of  F(Y), 
NS01A  initially  creates  a pseudo-cost  function  as 

J(Y)  = FT(Y)  F(Y)  (4-5) 

where  "T"  denotes  the  transpose  operation.  This  expression 
represents  a vector  multiplication  of  a (1  x n)  vector  times 
an  (n  x 1)  vector  which  produces  a quadratic  scalar  cost  J» 
where  J > 0 at  all  values  of  Y. 

A combination  of  either  Newton-Raphson  or  gradient  type 
iteration  steps  are  employed  to  find  the  minimum  of  the  func- 
tion J(Y)  to  within  a certain  accuracy  (ACC),  selected  by  the 
user.  To  start  the  numerical  solution  process,  the  user  must 
supply  NS01A  with  an  initial  guess  for  Y.  Iteration  then 
proceeds  between  NS01A,  which  checks  for 

J(Y)  < ACC  (4-6) 

and  CALFUN,  which  computes  the  values  of  F(Y)  for  each  iter- 
ation. Since  the  minimum  of  J is  zero,  the  values  of  Y which 
minimize  (4-5)  are  the  solution  values  to  the  nonlinear  equa- 
tions (4-2). 

The  NS01A  algorithm  is  quite  efficient  in  that  it  em- 
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ploys  a Newton-Raphson  technique,  but  automatically  switches 
to  the  method  of  steepest  gradient  when  it  detects  possible 
divergence  of  the  Newton-Raphson  steps.  So  NS01A  essential- 
ly uses  a gradient  step  when  the  initial  guess  of  Y is  far 
from  the  actual  solution,  then  provides  quadratic  conver- 
gence via  Newton-Raphson  steps  when  nearer  the  solution. 

Several  difficulties  are  inherent  in  the  application  of 
a general  purpose  subroutine  to  this  specific  problem.  First, 
the  values  of  the  ¥ variables  must  be  about  the  same  order  of 
magnitude  in  order  to  ensure  convergence . Since  four  of  the 
variables  are  angles  expressed  in  radians,  and  is  a TOF 
expressed  in  seconds,  the  usual  value  of  Y^  is  many  orders 
of  magnitude  different  from  the  other  components  of  Y.  This 
necessitated  scaling  the  TOF  variable.  The  proper  choice  of 
the  scaling  factor  greatly  aided  convergence. 

The  most  critical  problem,  however,  in  the  application 
of  NS01A  to  the  trajectory  targeting  process,  was  its  sen- 
sitivity to  the  initial  guess  of  Y.  Unless  this  guess  was 
fairly  close  to  the  actual  solution,  convergence  would  not 
occur . 

The  problem  of  good  initial  guesses  was  overcome  by 
first  targeting  the  transfer  using  the  assumption  that  the 
bums  were  impulsive . This  method  generated  values  of  Y 
which  were  sufficiently  close  to  the  actual  values  necessary 
for  the  real  finite  bum  case,  thus  facilitating  convergence 
by  NS01A. 

Since  targeting  is  done  to  satisfy  the  actual  finite 
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thrusting,  a numerical  scheme  must  be  used  to  integrate  the 
resulting  nonlinear  vector  differential  equations  of  motion 
(Appendix  C,  Equation  C-7).  Here  again  a general  purpose 
subroutine  combination  called  SET/STEP  is  employed  for  this 
purpose.  Thus,  CALFUN  calls  upon  SET/STEP  to  integrate  the 
equations  of  motion  resulting  from  the  current  values  of  Y. 

The  user  must  supply  SET/STEP  with  the  integration  step- 
size.  Then  SET  initializes  the  differential  equations,  and 
STEP  integrates  them  one  step  at  a time.  A classical  Runge- 
Kutta  method  is  used  for  the  first  three  steps,  and  then  a 
fourth  order  Adams-Bashforth-Adams-Moulton  predictor  correc- 
tor scheme  is  applied  to  succeeding  points. 

NS01A  gave  convergence  in  as  few  as  20  iterations,  or 
at  the  most  90  iterations  depending  on  how  closely  the  im- 
pulsively targeted  values  of  Y were  to  the  actual  values  for 
the  particular  transfer  trajectory  being  targeted. 

Maximum  Payload  Missions 

The  hit  conditions  employed  in  CALFUN  are  not  unique. 
Indeed,  many  different  forms  were  tried,  and  their  conver- 
gence properties  compared,  before  selecting  the  set  used  in 
(4-3).  There,  constraining  the  inclination  of  the  transfer 
trajectory  was  simply  a convenient  method  of  sampling  the 
range  of  trajectories  possible  for  excess  energy  missions. 

When  targeting  a transfer  for  near  maximum  payload  (for 
the  given  AV^,  AVg),  constraining  the  amount  of  first  burn 
plane  change  is  too  restrictive.  The  range  of  possible  tra- 
jectories narrows  considerably  for  payloads  near  maximum; 
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and  reduces  to  one  at  some  limiting  value.  In  these  cases 
it  is  difficult  to  predict  what  iT  should  be. 

An  alternate  set  of  hit  conditions  is  used  for  target- 
ing missions  where  the  payload  is  near  or  at  maximum.  This 
set  also  has  good  convergence  properties,  and  does  not  con- 
strain iTt 


pi  ■ - Vx 

p2  - V7>  * vy 

p3  * (4‘7) 

F4  = e2(y)  - 0 

p5  - i2<p>  - i2 


The  first  three  hit  conditions  now  become  the  components 
of  the  velocity  vector,  These  component  values  are  cal- 

culated at  each  call  of  CALFUN  to  correspond  to  the  orienta- 
tion of  r2(Y),  generated  by  that  iteration  of  NS01A. 

This  set  of  hit  conditions  will  also  converge  equally 
well  for  lower  payloads,  but  the  range  of  possible  trajec- 
tories present  in  this  case  (by  not  specifying  tjjg)  will  pro- 
duce a non-unique  solution. 

Targeting  S-fee  Tm 

This  step  may  be  thought  of  as  selecting  the  IUS  launch 
time  for  proper  phasing  between  the  orbits.  Step  one  of  tar- 
geting solves  for  five  of  the  variables,  referenced  to  the 
local  frame.  Once  values  for  these  five  variables  are  fixed, 
the  last  two  constraints  are  then  functions  of  the  only  var- 
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iable  left  free  in  step  one  (the  launch  time),  leaving 


n2*tbl> 

(4-8) 

T*<V> 

(4-9) 

Thus,  two  independent  equations  in  one  unknown  now  re- 
main. Since  the  equations  are  independent,  a value  for  ttl 
cannot,  in  general,  be  found  which  will  satisfy  both  rela- 
tionships. At  this  point  then,  either  (4-8)  or  (4-9)  is 
satisfied  depending  on  the  type  transfer  involved. 

In  the  case  of  a three-dimensional  transfer,  satisfy- 

# 

ing  TA  without  also  satisfying  Cl2  would  be  meaningless.  So 
if  the  transfer  is  between  non-conlanar  orbits,  the  target- 
ing in  this  step  automatically  satisfies  (4-8)  and  leaves 
(4-9)  free.  If  the  transfer  is  between  coplanar  orbits, 
then  a t^  which  satisfies  (4-9)  is  found  instead. 

Two  points  are  worthy  of  emphasis  here.  First,  even  if 
the  parking  orbit  and  mission  orbit  are  coplanar,  the  trans- 
fer orbit  itself  can  be  non-coplanar . An  out-of-plane  trans 
fer  trajectory  is  often  possible  because  of  excess  energy. 

In  this  case,  a component  of  AV1  is  depleted  by  thrusting 
out  of  plane  (creating  a non-coplanar  transfer  orbit),  and 
a component  of  AV2  is  likewise  utilized  to  regain  the  ini- 
tial plane.  For  this  transfer,  (4-9)  is  still  satisfied. 

The  second  point  is  that  even  though  a coplanar  trans- 
fer  is  necessary,  rendezvous  (i.e.,  satisfaction  of  TA  ) 
may  not  be  required.  In  this  case,  a t^  must  still  be 
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selected  so  that  a reference  location  is  specified  for  final 
computation  of  corresponding  thrust  direction  angles.  If 
rendezvous  is  not  a necessary  part  of  the  mission,  this  may 
be  accomplished  by  simply  specifying  some  arbitrary  value 
for  £02 j such  as  iQ2  = 0°.  This,  in  turn,  generates  an  ar- 
bitrary TA(t^),  which  when  solved  yields  a t^.  Here  any 
value  of  t^  will  do  just  as  well  if  rendezvous  is  unneces- 
sary, and  it  gives  a reference  point  in  the  parking  orbit 
needed  to  compute  the  actual  thrust  direction  angles  ex- 
pressed in  the  geocentric-equatorial  frame. 

The  value  of  ttl,  once  determined,  represents  the  ab- 
solute time  in  seconds  after  the  epoch  time  when  the  first 
opportunity  for  transfer  occurs.  The  actual  time  for  second 
stage  ignition,  then,  now  becomes 

- *bl  ♦ WF*.  tbb  (“-10) 

where  T0F*is  the  total  time  in  seconds  from  first  stage  ig- 
nition to  second  stage  burnout  (as  previously  computed  in 
targeting  step  one),  and  t^b  is  the  fixed  burn  time  of  the 
second  stage. 

The  next  operation  carried  out  in  step  two  of  targeting 
is  to  translate  the  values  of  the  thrust  direction  angles 
expressed  in  the  local  frame  from  targeting  step  one,  to 
their  values  expressed  in  the  geocentric-equatorial  frame. 
This  can  be  done  through  appropriate  coordinate  transforma- 
tions once  t^  is  known.  As  explained  earlier,  given  the 
orbital  elements  of  the  parking  orbit,  ric(t)  may' be  computed 
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in  the  geocentric-equatorial  frame.  When  is  known,  the 
inertial  position  for  first  stage  ignition  is  fixed  as  ric 
(t^i).  That  vector,  in  turn,  fixes  the  orientation  between 
the  local  frame  and  the  geocentric-equatorial  frame,  as 
shown  in  Figure  4-4,  for  the  example  of  a geosynchronous 
transfer  mission. 

The  last  step  in  targeting  i:s  to  compute  the  next  four 
sequential  times  when  the  mission  transfer  could  occur,  and 
the  corresponding  thrust  direction  angles  for  each  time.  The 
reason  for  this  is  so  these  values  could  be  stored  on-board 
the  IUS  and  sequentially  used  as  necessary  if  mission  delays 
occur.  Here  the  usefulness  of  doing  the  original  targeting 
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in  a local  frame  is  apparent.  When  each  of  the  next  four 
tfel  times  is  determined,  it  is  only  necessary  to  go  through 
the  same  coordinate  transformation  for  each  new  ric^bl^ 
position,  rather  than  retarget  the  whole  mission  again.  The 
fact  that  five  sequential  times  are  chosen  here  is  completely 
arbitrary.  Any  number  of  parameter  value  sets  could  equally 
well  be  calculated  in  this  way. 

If  the  transfer  is  between  non-coplanar  orbits,  the 
interval  between  successive  mission  opportunity  times  is 
just  equal  to  the  period  of  the  parking  orbit,  and  the  thrust 
direction  angles  remain  the  same.  This  is  because  the  trans- 
fer must  be  initiated  at  only  one  position  in  the  transfer 
orbit.  If  the  transfer  is  between  coplanar  orbits,  then 
the  interval  between  successive  ttl  times  is  just  equal  to 
the  synodic  period  of  the  two  orbits,  and  the  values  of  the 
thrust  direction  angles  are  different  at  each  opportunity. 

The  actual  calculations  discussed  in  this  section  may 
be  found  in  Appendix  D. 


where  AM  is  a (6  x l)  "miss"  vector  (insertion  error  vector), 
AAL  is  a (3  x l)  thrust  misalignment  vector,  and  (BM/dAL)  is 
a (6  x 6)  matrix  of  sensitivity  coefficients.  Written  out 
in  full  this  expression  is 
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where  the  subscript  "i"  denotes  "at  insertion".  All  the  in- 
dividual components  and  sensitivities  are  expressed  with  re- 
spect to  the  local  frame.  This  enables  a more  meaningful 
comparison  of  results  between  diverse  transfers.  The  units 
of  the  individual  sensitivities  are  nm  of  position  error  per 
milliradian  of  thrust  misalignment!  or,  ft/sec  of  velocity 
error  per  milliradian  of  thrust  misalignment. 

With  the  sensitivities  in  this  form,  however,  it  is 
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difficult  to  draw  a standard  of  comparison  for  evaluating 
the  accuracies  of  individual  trajectories!  to  determine 
which  one  from  a group  of  possible  trajectories  would  give 
the  lowest  insertion  errors.  A great  deal  of  thought  went 
into  how  best  to  use  the  information  provided  by  the  thrust 
misalignment  sensitivity  matrices  to  compare  trajectory 
accuracies. 

This  dilemma  was  finally  resolved  by  deciding  upon  a 
"worst  case"  comparison.  But  before  explaining  this  method 
it  is  instructive  to  look  at  the  dilemma  in  the  light  of 
what  the  sensitivities  mean.  The  overall  miss  vector  (AM) 
can  be  broken  into  two  components;  one  for  position  error, 
Ar,  and  the  second  for  velocity  error,  AV,  as  the  partition- 
ing in  (5-2)  shows.  Taking  for  example  the  position  error 
vector  Ar  (same  for  velocity),  the  error  is  expressible  as 


Ar  = Arxii  + Aryij  + Arz-k 


/drxi 
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(5-3) 


where  the  magnitude  of  the  position  error  is 


Ar  = 7a  rxi2  + Aryi2  + Arzi2 

Equation  (5-4)  in  turn  depends  on  the  actual  components 


(5-4) 
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of  the  misalignment  vector 


AAL  = AL  i + AL  j + AL_  k ( 5~5) 

x y l 

Thus,  as  seen  from  studying  Equation  (5-3) » the  actual  value 
of  Ar  for  any  particular  transfer  is  dependent  on  the  par- 
ticular magnitudes  (and  even  the  signs)  of  the  individual 
components  of  AAL.  But  AAL  is  a completely  random  variable! 
At  best,  only  a tolerance  (or  upper  limit)  on  the  magnitude 
of  AAL  may  be  known. 

Assuming  an  arbitrary  orientation  and  magnitude  for  AAL, 
and  using  this  same  value  to  compute  a subsequent  Ar  for  each 
trajectory,  is  one  possible  approach.  However,  this  approach 
is  still  inadequate  in  that  there  is  no  way  of  knowing  wheth- 
er or  not  this  arbitrary  AAL  has  the  same  effect  on  each  of 
the  widely  diverse  trajectories.  Also,  there  would  be  no 
way  of  predicting  whether  a "worse”  AAL  direction  might  be 
possible.  This  approach,  however,  is  meaningful  if  a guar- 
anteed worst  alignment  direction  were  used  based  on  the  sen- 
sitivity characteristics  of  each  trajectory. 

Although  there  is  no  way  of  predicting  the  actual  orien- 
tation of  the  random  thrust  misalignment  vector,  it  is  valid 
to  say  that  there  will  be  in  every  case  a "worst  possible" 
orientation.  That  is,  there  will  be  some  orientation  of  the 
AAL  vector,  which  (for  any  fixed  magnitude)  will  cause  the 
greatest  insertion  errors.  This  concept,  then,  is  imple- 
mented here. 

The  sensitivity  matrix  of  each  transfer  is  reduced  to 
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two  "worst  case”  sensitivities;  one  for  position  error  and 
one  for  velocity  error.  To  explain  how  the  worst  case  sen- 
sitivities are  found,  the  example  of  insertion  position  er- 
ror can  again  be  used.  Derivation  of  the  worst  case  veloci- 
ty error  sensitivity  is  exactly  the  same.  Using  a more  sim- 
plified notation  to  convey  the  concept,  let 

y = A x (5-6) 

where  y = Sr,  the  insertion  position  error  vector;  x = AAL, 
the  thrust  misalignment  vector;  and  A is  the  matrix  of  sen- 
sitivity coefficients,  (dr-/dAAL),  which  transforms  the  vec- 
tor x into  the  vector  y. 

Now  the  question  of  a worst  misalignment  direction  may 
be  put  in  this  form;  For  a fixed  specified  magnitude  of  the 
vector  x,  what  is  the  maximum  possible  magnitude  of  the  vec- 
tor y.  Since  a sensitivity  is  desired,  this  equates  to  ask- 
ing what  orientation  of  the  unit  thrust  misalignment  vector 
goes  through  the  vector  transformation  of  (5-6)  to  give  the 
maximum  length  of  the  insertion  position  error  vector. 

With  the  assumption  that  the  thrust  misalignment  is 

strictly  due  to  misalignment  between  the  XYZ  frame  and  the 

X Y Z frame,  the  worst  possible  orientation  of  the  thrust 
P P P 

misalignment  vector  is  along  the  "Euler  angle"  axis  direc- 
tion which  produces  the  "worst"  misalignment  between  these 
two  frames*  in  the  sense  of  causing  maximum  insertion  error. 
This  is  illustrated  in  Figure  5-1. 

As  derived  in  Appendix  E,  the  maximum  possible  magni- 
tude of  y caused  by  some  x is  obtained  from 
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where  k is  the  largest  eigenvalue  of  the  A A matrix. 

Thus,  by  computing  the  eigenvalues  of  the  symmetric 
matrix  ATA  formed  from  the  respective  position  or  velocity 
error  sensitivity  matrices,  the  two  worst  case  sensitivi- 
ties are  obtained.  From  these,  and  a value  for  the  magni- 
tude of  the  worst  thrust  misalignment  expected,  the  upper 
bounds  of  the  associated  insertion  errors  are  obtained  di- 
rectly. These  values  may  then  be  compared  between  an  assort- 
ment of  trajectories  as  a valid  and  direct  measure  of  their 
respective  accuracies. 

Thrust  Magnitude  Error  Sensitivities 

Vehicle  performance  is  calculated  using  the  ideal  ve- 
locity UV)  equation 
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(2-1) 


where  I is  the  specific  impulse  of  the  solid  propellant; 
gQ  is  the  gravitational  constant  (32.14644  ft/sec  ) ; mQ  is 
the  vehicle  mass  at  ignition;  and  mf  is  the  vehicle  mass  at 
burnout.  From  (2-1)  it  can  be  seen  that  the  AW  provided  by 
a specific  stage  is  a function  only  of  the  propulsion  pa- 
rameters (I__  and  propellant  weight),  and  the  inert  (empty) 

sp 

weight  of  the  vehicle. 

It  is  assumed  here  that  fairly  tight  tolerances  can  be 

maintained  on  I , fuel  weight  and  the  structure  weight  of 
sp 

the  vehicle  so  that  the  magnitude  of  AW  for  each  stage  is 
not  subject  to  error.  However,  a deviation  in  the  thrust 
profile  of  either  engine  is  possible  without  any  variation 
in  its  total  impulse  ( AW ).  That  is,  the  total  velocity 
change  can  be  equal  to  that  expected,  but  the  acceleration 
profile  may  vary.  This  would  be  caused  by  a variation  in 
the  solid  propellant  bum  rate,  thus  causing  a deviation  in 
the  thrust  magnitude.  This  can  be  seen  from  the  thrust 
equation 


T = -cm 


(5-8) 


where  c is  the  effective  exhaust  velocity;  and  m is  the  mass 
flow  rate.  The  mass  flow  rate  equation  (assuming  constant  m) 


m(t)  = mQ  - mt 


(5-9) 


Evaluating  this  equation  at  the  burnout  time  (t^),  and  solv- 
ing for  t^,  gives  the  total  bum  time  as 


m0  ~ mf 
^b  “ in 


(5-10) 


So  a variation  in  the  magnitude  of  T will  be  associated  with 
a change  in  the  bum  time.  A variation  in  the  bum  time 
will  affect  the  position  at  burnout.  This  can  be  a major 
source  of  insertion  error  since  (as  discussed  in  Chapter  I) 
the  required  velocity  to  hit  the  target  is  a function  of  the 
burnout  position. 

As  with  thrust  direction  error,  the  thrust  magnitude 
error  is  in  general  a random  variable.  However,  bounds 
on  this  error  should  be  predictable  from  studying  the  sta- 
tistical characteristics  of  the  engine  performance.  With 
this  in  mind,  sensitivities  are  calculated  for  both  an  over- 
all plus  and  an  overall  minus  lf>  error  in  the  thrust  magni- 
tude. For  each  sensitivity  it  is  assumed  that  both  engines 
have  the  same  thrust  error. 


VI.  Results  for  the  Reference  Missions 


This  chapter  presents  the  results  of  targeting  and  the 
accuracy  analysis  of  the  geosynchronous  and  subsynchronous 
reference  missions.  The  geosynchronous  mission  is  given  the 
greatest  emphasis.  The  results  for  this  mission  are  present- 
ed for  two  different  fuel  load  combinations}  each  with  four 
different  payload  weights.  In  addition,  the  results  for  a 
coplanar  (geosynchronous)  transfer  are  included  - along  with 
an  expanded  explanation  of  some  of  the  intermediate  targeting 
and  error  analysis  steps  - since  this  maneuver  is  easily  vis- 
ualized. 

For  targeting  purposes,  arbitrary  values  were  assumed 
for  Z02  (true  longitude  at  epoch  of  the  mission  orbit  target 
position),  and  0-^  (longitude  of  the  ascending  node  of  the 
parking  orbit).  The  values  obtained  for  the  six  control  pa- 
rameters for  each  mission  are  not  included  here,  since  they 
are  based  on  arbitrary  initial  conditions  and  the  selected 
plane  change  split  between  bums. 

The  results  of  the  error  analyses  are  presented  in  graph 
ical  form  for  most  missions.  Insertion  error  sensitivities 
are  plotted  against  the  first  bum  plane  change  angle,  for 
the  various  fuel  and  payload  combinations.  The  sensitivities 
due  to  thrust  misalignment  are  the  "worst  case"  values.  The 
units  on  the  position  and  velocity  insertion  error  sensitivi- 
ties due  to  thrust  misalignment  are  natuical  miles  and  ft/sec 
per  milliradian  of  thrust  vector  misalignment;  and  similarly 
for  insertion  error  sensitivities  due  to  thrust  magnitude 
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variations,  per  ± 1 % variation.  The  last  plot  for  each  mis- 
sion shows  the  transfer  time  of  flight  versus  first  bum  plane 
change  angle. 

In  all  cases  except  two,  the  maximum  value  of  the  first 
bum  plane  change  angle  shown  on  the  plots  is  within  1°  of 
its  absolute  maximum.  The  two  exceptions  are  shown  in  Fig- 
ures 6-7,  6-8,  6-9.  and  in  Figures  6-16,  6-17,  6-18,  where 
the  maximum  value  of  the  first  bum  plane  change  angle  was 
iterated  to  within  .1°  of  its  absolute  maximum. 

The  last  section  in  this  chapter  summarizes  the  overall 
results,  and  gives  some  general  observations  concerning  this 
scheme . 

IUS  Vehicle  Specifications 

Table  III  lists  the  values  of  the  IUS  vehicle  parameters 
which  were  used  in  targeting  the  reference  missions.  These 
values  are  based  (with  some  minor  simplifications)  on  the 
Burner  II  specifications  as  given  in  Reference  2. 


Table  III.  Vehicle  Specifications 


First  Stages 

Total  Structure  Weight  (ST^ 

2437  lb 

Propellant  Weight  (PR0P1) 

17,300  lb  or 
20,000  lb 

Thrust  (T^) 

41,923.4  lbf 

Specific  Impulse  (Ispl) 

291.8  sec 

Second  Stage 

Total  Structure  Weight  (ST2) 

1362  lb 

Propellant  Weight  (PR0P2) 

4700  lb 

Thrust  (T2) 

14,345.6  lbf 

Specific  Impulse  (ISp2) 

300.8  sec 
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The  Geosynchronous  Mission 

Mission  Description:  the  parking  orbit  is  circular,  at 

an  altitude  of  160  nm  and  is  inclined  28. 5°.  The  mission 
orbit  is  equatorial,  and  at  an  altitude  of  19,323  nm. 

Figures  6-1  through  6-9  show  the  results  for  a fuel  load 
combination  of  PROP^  = 17,300  lbs,  PR0P2  = 4700  lbsj  and  pay- 
load  (PL)  weights  of  1000,  2000,  and  3000  lbs.  Similarly, 
Figures  6-10  through  6-18  show  the  comparable  results  for 
a fuel  loading  of  PROP-,^  = 20,000  lbs,  PR0P2  = 4700  lbs. 

To  obtain  results  for  payloads  very  near  maximum,  the 
alternate  hit  conditions  given  in  Equation  4-7  were  employed. 
In  these  cases,  the  amount  of  first  bum  plane  change  is  left 
free  to  seek  its  optimum  value.  Table  IV  lists  the  maximum 
payload  results  for  the  two  different  fuel  load  combinations. 

The  retargeting  period  is  5427  seconds  for  all  the  geo- 
synchronous non-coplanar  transfer  missions. 

The  Subsynchronous  Mission 

Mission  Description:  The  parking  orbit  is  at  an  alti- 

tude of  160  nm,  and  inclined  57° • The  mission  orbit  is  at 
an  altitude  of  10,900  nm,  and  has  an  inclination  of  63°. 

The  longitudes  of  the  ascending  nodes,  and  02,  are  equal. 

The  subsynchronous  transfer  mission  will  be  used  for 
placement  of  the  Department  of  Defense  Navstar  Global  Posi- 
tioning System  (GPS).  The  total  payload  weight  is  forecast 
to  be  4400  lbs.  The  results  for  a subsynchronous  mission 
with  this  amount  of  payload  are  shown  in  Figures  6-19,  6-20, 
and  6-21.  The  fuel  loading  for  this  transfer  was  PR0P1  = 


PL  = 1000  LBS 


Thrust  Magnitude 


Thrust  Direction 


2 4 6 8 10 

FIRST  BURN  PLANE  CHANGE  (DEG) 


Figure  *-l.  Geosynchronous  Insertion  Error 
Sensitivities  (Position) 

(PROP1  = 17.300  lb,  PL  = 1000  lb) 


INSERTION  VELOCITY  ERROR  IFT/SEC) 


E OF  FLIGHT  I SEC) 


INSERTION  POSITION  ERROR  (NM) 


Figure  6-4.  Geosynchronous  Insertion  Error 
Sensitivities  (Position) 

(PROPj^  = 17,300  lb,  PL  = 2000  lb) 
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40-1 


PL  = 2000  LBS 


Figure  6-5.  Geosynchronous  Insertion  Error 
Sensitivities  (Velocity) 

(PR0P1  = 17,300  lb,  PL  = 2000  lb) 
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PL  = 3000  LBS 
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Figure  6-7. 


Geosynchronous  Insertion  Error 
Sensitivities  (Position) 

(PROPj^  = 17,300  lb,  PL  = 3000  lb) 
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PL  = 3000  LBS 
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Figure  6-8.  Geosynchronous  Insertion  Error 
Sensitivities  (Velocity) 

(PROP1  = 17,300  lb,  PL  = 3000  lb) 


69 


Figure  6-9.  Geosynchronous  Time  of  Flight 

(PROP1  = 17,300  lb,  PL  = 3000  lb) 
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PL  = 1000  LBS 
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Figure  6-10.  Geosynchronous  Insertion  Error 
Sensitivities  (Position) 

PROP,  = 20,000  lb,  PL  = 1000  lb 


15000-1 


PL  = 1000  LBS 


Figure  6-12.  Geosynchronous  Time  of  Flight 

(PROP-l  = 20,000  lb,  PL  = 1000  lb) 
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Figure  6-13.  Geosynchronous  Insertion  Error 
Sensitivities  (Position) 

( PROP x = 20,000  lb,  PL  » 2000  lb) 
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TRANSFER  TIME  OF  FLIGHT  (SEC) 
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PL  = 2000  LBS 
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Figure  6-15.  Geosynchronous  Time  of  Flight 

(PROP-j^  * 20,000  lb,  PL  = 2000  lb) 
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Figure  6-1 6.  Geosynchronous  Insertion  Error 
Sensitivities  (Position) 

(PROP,  = 20,000.1b,  PL  = 3000  lb) 
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PL  = 3000  LBS 
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Figure  6-17. 


Geosynchronous  Insertion  Error 
Sensitivities  (Velocity) 


(PROP1  = 


20,000  lb,  PL  = 3000  lb) 
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Figure  6-18,  Geosynchronous  Time  of  Flight 

(PROP-,  = 20,000  lb,  PL  = 3000  lb) 


PL  = 4400  LBS 
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Figure  6-20, 


Subsynchronous  Insertion  Error 
Sensitivities  (Velocity) 

(PR0P1  * 17,300  lb,  PL  = **400  lb) 


TRANSFER  TIME  OF  FLIGHT  l SEC) 


Figure  6-21.  Subsynchronous  Time  of  Flight 

(PROP1  = 17,300  lb,  PL  = 4400  lb) 
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Table  IV.  Maximum  Payload  Geosynchronous  Transfers 


Fuel  Loading  (lbs) 

PROP 2=17, 300 
PROP 2=4700 

PR0P1=20, 000 
PROP 2=4700 

Payload 

4250 

4600 

First  Bum  Plane 
Change  (deg.) 

1.7 

5.9 

Insertion  Error. 
Position  (nm/mr) 

11.4 

36.0 

Insertion  Error, 
Velocity  (ft/sec/mr) 

9.1 

21.5 

Position  Ins.  Error, 
± 1%  Thrust  Dev. 

21.0 

32.4 

Velocity  Ins.  Error 
± ljt  Thrust  Dev. 

10.5 

17.3 

Time  of  Flight 
(sec) 

18,954 

18,440 

Transfer  Angle 
(deg. ) 

179.9 

175.2 

17,300  lbs,  PR0P2  = 4700  lbs. 

The  retargeting  period  for  this  mission  is  5427  seconds, 
the  same  as  for  the  geosynchronous  mission. 

Transfer  Between  Coplanar  Orbits 

For  a transfer  between  any  two  (circular)  coplanar  or- 

bitB,  an  in-plane  transfer  trajectory  can  always  be  found 

which  matches  any  allowable  and  aV2  combination. 

With  an  IUS  vehicle  propellant  and  payload  combination 

of  PROPj^  = 20,000  lbs,  PR0P2  = 4700  lbs  and  PL  = 3000  lbs, 

the  consequent  AV  capabilities  become  9443  ft/sec  for  AV, , 
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and  7070  ft/sec  for  aV2.  As  shown  in  Figure  3-4,  this  is  an 
allowable  combination  to  complete  a coplanar  transfer  from  a 
160  nm  parking  orbit  to  a geosynchronous  mission  orbit.  Be- 
cause a coplanar  transfer  is  relatively  easy  to  visualize,  the 
applicable  energy  management  (non-Hohmann)  trajectory  for  this 
transfer  will  be  briefly  described  here. 

The  impulsive  trajectory  which  matches  this  AV  combina- 
tion is  one  with  a transfer  angle  (central  angle)  of  140°,  a 
first  bum  flight  path  angle  of  8.6°  and  a second  bum 

flight  path  angle  (4^)  of  44.5°.  Referenced  to  the  local 
inertial  frame  (as  the  transfer  appears  in  Figure  3-1).  im- 
pulsive targeting  yields  a first  bum  thrust  direction  angle 
(t^)  of  57.7°,  and  a second  bum  thrust  direction  angle  (p-j) 
of  -82.9°.  The  impulsive  TOF  is  12,107  seconds. 

When  the  targeting  is  refined  to  include  the  finite  bum 
dynamics,  the  transfer  angle  becomes  145.2°,  and  the  thrust 
direction  angles  become  ^ = 62.8°,  = -77. 8° 1 with  a 

transfer  TOF  of  12,230  seconds.  It  is  interesting  to  note 
that  in  this  case,  as  well  as  nearly  all  the  cases  targeted, 
the  differences  in  the  thrust  angles  between  their  values 
using  the  impulsive  approximation  and  their  actual  finite 
bum  values  is  about  5°*  This  is  a significant  difference, 
and  demonstrates  that  open  loop  targeting,  without  including 
finite  bum  dynamics,  would  be  wholly  inadequate. 

The  accuracy  analysis  of  this  trajectory  produced  a 
position  and  velocity  insertion  error  sensitivity  matrix 
(dM/aKE  matrix  of  Equation  5-2),  due  to  thrust  misalignment, 
with  the  following  values 1 
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r 


aM 

.0047 

-.0195 

3.6696 

-.0019 

-.0019 

-1.9258 

-12.1966 

12.6168 

0.0 

a At 

.0080 

-.0010 

-2.4489 

-.0092 

- . 0004 

12.7996 

-7 . 0846 

-1.3990 

0.0 

where  the  units  are  nm/mr,  and  (ft/sec)/mr . Combining  these 
into  "worst  case"  sensitivities  via  Equation  5-7  gave  an  in- 
sertion position  error  of  17.548  nm/mr,  and  an  insertion  ve- 
locity error  of  13.032  ( ft/sec )/mr. 

The  insertion  error  sensitivities  due  to  thrust  magni- 
tude deviations  are  31.331  run  and  20.936  ft/sec  for  positive 
deviations » and  31.974  nm  and  21.382  ft/sec  for  negative  devi- 
ations. 

The  synodic  period  (interval  between  transfer  opportuni- 
ties to  effect  a rendezvous)  for  this  coplanar  transfer  is 
5791  seconds. 

Observations  Concerning  the  Results 

Overall,  the  results  obtained  by  this  scheme  are  very 
encouraging.  It  is  clear  that  for  most  energy  management 
transfers  a significant  amount  of  targeting  flexibility  is 
available  to  the  mission  planner,  allowing  an  "optimal"  tra- 
jectory to  be  selected  from  a wide  range  of  possibilities. 

The  sensitivities  to  error  inputs  and  transfer  times  vary 
significantly  over  the  range  of  usable  transfers.  Thus,  the 
final  choice  of  a trajectory  would  be  somewhat  vehicle  and 
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mission  dependent. 

The  plots  show  that  insertion  errors  caused  by  thrust 
misalignment  always  increase  as  the  amount  of  plane  change 
accomplished  by  the  first  bum  (i|j2)  increases.  However,  in- 
sertion errors  due  to  thrust  magnitude  deviations  always  de- 
crease with  an  increase  in  vj/g • As  a general  rule,  insertion 
errors  increase  with  an  increase  in  fuel  and/or  payload.  The 
transfer  time  of  flight  increases  with  ijjgi  with  the  minimum 
always  occurring  where  iJj2  = 0. 

In  order  to  be  meaningful,  insertion  error  sensitivities 
must  be  multiplied  by  the  expected  amount  of  error  input.  En- 
gineers currently  working  on  IUS  design  indicate  that  the  up- 
per limit  on  IMU  misalignment  prior  to  the  first  bum  is  1.8 
mr,  and  that  the  upper  bound  on  expected  thrust  magnitude  var- 
iations is  ± 

The  required  IUS  orbital  insertion  accuracies  are  given 
in  Reference  5.  That  document  specifies  a maximum  permissi- 
ble position  error  of  ± 92  nm,  and  maximum  allowable  velocity 
error  of  ± 78  ft/sec. 

Assuming  the  only  error  inputs  to  be  those  caused  by  er- 
rors in  the  application  of  the  thrust  vectors,  trajectories 
can  be  selected  which  satisfy  the  insertion  error  requirements 
for  each  of  the  missions  targeted.  Por  example,  a geosynch- 
ronous transfer  with  17,300  lbs  of  first  stage  propellant  and 
3000  lbs  of  payload  is  a good  baseline  for  comparison.  Plots 
for  this  mission  are  shown  in  Figures  6-7,  6-8  and  6-9.  By 
using  a transfer  trajectory  with  0°  of  first  bum  plane  change, 
and  asstuning  the  most  pessimistic  case  in  which  insertion 
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errors  due  to  thrust  misalignment  and  thrust  magnitude  var- 
iations would  add  linearly,  the  resulting  insertion  errors 
would  be  36  nm  and  26  ft/sec.  These  figures  are  worst  case 
values  based  on  the  upper  bound  values  for  the  error  inputs 
mentioned  earlier.  These  errors  are  well  within  the  allow- 
able tolerances  of  ± 92  nm  and  ±78  ft/sec.  In  addition,  this 
trajectory  (with  = 0°)  yields  the  minimum  transfer  time. 

A minimum  time  of  flight  trajectory  would  be  desirable  to 
minimize  the  IMU  platform  gyro  drift  prior  to  the  second  burn. 

It  should  be  emphasized  that  thrust  vector  errors  may 
not  be  the  only  significant  error  inputs  affecting  the  trans- 
fer trajectory.  Also  present  to  some  extent  will  be  inser- 
tion errors  due  to  Isp  and  vehicle  (fuel  and  inert)  weight 
dispersions.  These  dispersions  can  cause  variations  in  the 
magnitude  of  the  velocity  impulse  applied  by  each  bum  (ref- 
erence Equation  2-1),  which  was  assumed  invariant  for  this 
study. 

Another  important  consideration  in  interpreting  the  sen- 
sitivities is  determination  of  a realistic  multiplication 
factor  for  overall  thrust  misalignment,  which  would  take  in- 
to account  any  gyro  drift  between  bums.  If,  however,  a re- 
alignment of  the  IMU  platform  could  be  accomplished  using  a 
star  tracker  or  similar  means  prior  to  the  second  bum,  the 
question  of  gyro  drift  would  not  apply.  This  would  leave 
the  velocity  impulse  errors  (as  described  above)  as  the  only 
other  important  error  source  yet  to  be  applied  to  this  scheme 
(assuming  negligible  modeling  errors). 

Staying  with  the  concept  of  simplicity  in  design,  it  may 
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be  possible  to  accomplish  the  burns  open  loop  also,  as  op- 
posed to  closed  loop  (constant)  attitude  control  during  the 
bums.  This  would  negate  the  requirement  for  a Thrust  Vector 
Control  (TVC)  system,  and  any  associated  on-board  software. 

The  feasibility  of  open  loop  bums  would  depend  on  how  ac- 
curately a fixed  nozzle  design  could  keep  the  applied  thrust 
vector  aligned  with  the  vehicle  center  of  gravity,  as  any 
significant  misalignment  would  cause  some  vehicle  attitude 
rotation  during  the  bums.  Here  again,  a realistic  multi- 
plication factor  could  be  computed  which  would  include  the 
expected  thrust  vector  misalignment  caused  by  any  vehicle 
rotation  during  the  bums.  This  multiplication  factor  could 
be  used  in  conjunction  with  the  thrust  misalignment  sensi- 
tivities computed  here,  to  indicate  the  insertion  error  that 
would  be  associated  with  open  loop  bums. 

A final  point  of  interest  concerning  mission  flexibil- 
ity should  be  made.  The  range  of  transfers  available  for  a 
particular  mission  allows  a wide  range  of  possible  transfer 
angles,  which  adds  a limited  degree  of  freedom  to  non-co- 
planar  transfers  which  involve  a rendezvous  maneuver.  That 
is,  if  the  phasing  between  the  parking  orbit  and  mission  or- 
bit was  approximately  correct,  the  flexibility  in  transfer 
angles  available  would  allow  a transfer  to  be  chosen  which 
"matches"  the  exact  phasing  for  rendezvous. 

Using  the  geosynchronous  transfer  example  cited  earlier, 
which  involved  a first  stage  propellant  weight  of  17,300  lbs 
and  a payload  of  3000  lbs,  the  range  of  admissible  values 
for  the  first  bum  plane  change  angle  (^2)  is  from  0°  to  5.9°. 


88 


The  associated  range  of  transfer  angles  is  157.0  - 172.3  . 

With  a first  stage  fuel  load  of  20,000  lbs.  the  range  of 
increases  to  0°  - 9.9°.  and  the  transfer  angle  range  becomes 
148.9°  - 176. 3° > a 27.4°  spread. 

In  addition  to  the  missions  presented  in  this  chapter, 
many  other  possible  transfers  were  targeted  for  various  or- 
bits, vehicle  parameters,  and  payload  weights.  In  all  cases 
the  results  obtained  were  consistent  with  the  results  for  the 
two  reference  missions  shown  here. 
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VII.  Conclusions  and  Recommendations 

Trajectory  matching  appears  to  be  a highly  effective  en- 
ergy management  technique . By  selecting  a non-Hohmann  tra- 
jectory which  matches  the  fixed  energy  capabilities  of  the 
IUS  vehicle,  a wide  range  of  possible  transfer  trajectories 
is  available  for  any  mission.  A great  deal  of  flexibility  is, 
therefore,  available  to  the  mission  planner,  allowing  an  "op- 
timal" transfer  trajectory  to  be  selected. 

The  use  of  a nonlinear  equation  solving  routine  was  ex- 
tremely advantageous  in  that  it  allowed  targeting  to  be  ac- 
complished to  include  the  finite  burn  dynamics,  without  any 
inherent  algorithm-generated  insertion  error.  This  exact 
targeting  is  a necessary  part  of  the  open  loop  design. 

Using  upper  bounded  values  for  the  thrust  vector  error 
inputs,  the  accuracy  analysis  showed  that  the  insertion  er- 
rors were  in  all  cases  well  within  the  allowable  limits.  The 
most  significant  insertion  error  source  turned  out  to  be 
thrust  vector  misalignment.  For  this  case,  the  trajectory 
which  minimizes  the  insertion  errors  is  also  the  minimum  time 
of  flight  trajectory.  This  trajectory  is  obtained  by  a plane 
change  split  between  the  burns  which  places  the  minimum 
amount  of  plane  change  (usually  zero)  in  the  first  burn. 
Insertion  errors  caused  by  thrust  misalignment  always  in- 
crease with  increasing  first  burn  plane  change  angle,  where- 
as the  insertion  errors  caused  by  thrust  magnitude  deviations 
always  decrease  with  an  increase  in  the  first  bum  plane 
change  angle. 
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The  results  of  this  study  show  that  open  loop  guidance  is 
quite  feasible  if  the  additional  effects  of  error  inputs  other 
than  thrust  vector  perturbations  can  be  controlled.  In  this 
regard,  gyro  drift  errors  could  be  effectively  eliminated  by 
a star  tracker  realignment  prior  to  the  second  burn.  If  in- 
sertion errors  due  to  Igp  and  vehicle  (fuel  and  inert)  weight 
dispersions  are  found  to  be  significant,  it  may  be  more  ec- 
onomical in  the  long  run  to  attempt  to  achieve  tighter  tol- 
erances on  these  dispersions  and  utilize  the  simple  open  loop 
design,  rather  than  invest  in  closed  loop  software  and  RCS 
correction  bums. 

It  is  recommended  that  this  stt  .j  be  extended  to  include 
targeting  for  transfers  between  elliptical  orbits.  In  ad- 
dition, a more  rigorous  analysis  of  this  scheme  could  be  ac- 
complished by  the  addition  of  the  remaining  error  sources 
(primarily  gyro  drift  and  velocity  impulse  errors)  followed 
by  a Monte  Carlo  type  statistical  analysis  of  the  subsequent 
insertion  errors. 
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Appendix  A 

Computer  Simulation  Algorithm 


A verbal  flow  chart  is  given  here,  which  summarizes  the 
sequence  of  steps  performed  by  the  computer  program  used  for 
this  study. 

Targeting  Portion 

1.  From  the  input  data  (as  listed  in  Chapter  II),  com- 
putes available  AV^  and  AV2j  corrects  AV^  for  estimated  fi- 
nite bum  losses. 

2.  Computes  the  range  of  first  bum  plane  change  pos- 
sible . 

3.  Targeting  step  one:  Sets  the  amount  of  first  bum 

plane  change  to  be  accomplished  to  the  lower  limit,  and  us- 
ing the  impulsive  approximation,  targets  a transfer  trajec- 
tory to  match  AV^  and  AV2,  which  satisfies  constraints  i<p , 

r2,  V2,  i2  and  e2>  Results  are  values  for  thrust  direction 

# * 

angles,  transfer  TOF  , and  transfer  angle  TA  , all  expressed 
in  the  local  frame. 

4.  Using  the  nonlinear  equation  solving  subroutine 
(NS01A) . with  its  initial  guesses  as  the  values  from  the  im- 
pulsive targeting  of  the  last  step,  and  numerical  integration 
of  the  equations  of  motion,  solves  for  the  actual  finite  bum 
values  of  the  control  parameters,  still  referenced  to  the 
local  frame. 

5.  Targeting  step  twoi  From  the  specified  orbital  el- 

« 

ments  of  each  orbit,  computes  ric(t)  and  r2C(t+T0F  ),  and 
from  these  determines  (first  mission  opportunity  time 
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after  epoch),  to  satisfy  either  02  if  the  transfer  is  between 

* 

non-coplanar  orbits,  or  to  satisfy  TA  (rendezvous)  if  the 
transfer  is  between  coplanar  orbits. 

6.  Transforms  the  values  for  the  thrust  direction  an- 
gles (expressed  in  the  local  frame)  to  their  corresponding 
values  in  the  geocentric-equatorial  frame.  Outputs  are  the 
actual  values  for  the  six  control  parameters  which  would  be 
stored  on-board  the  IUS. 

?.  Mission  Delay  retargeting!  Computes  the  values  of 
the  six  targeting  constants  which  correspond  to  the  next  four 
sequential  mission  opportunity  times,  using  appropriate  co- 
ordinate transformations . 


8.  Computes  the  insertion  position  and  velocity  error 
sensitivity  coefficient  matrices  due  to  thrust  misalignment, 
for  the  target  trajectory. 

9.  Reduces  these  matrices  to  two  "worst  case"  sensi- 
tivities, one  for  insertion  position  error  and  one  for  in- 
sertion velocity  error. 

10.  Computes  insertion  position  and  velocity  error  sen- 
sitivities due  to  thrust  magnitude  fluxuations. 

11.  Increases  the  value  of  the  first  burn  plane  change 
by  1°,  then  returns  to  (1)  and  repeats  steps  (1-10).  Contin- 
ues iterating  until  all  the  transfer  trajectories  in  the 
range  of  possible  first  burn  plane  change  have  been  targeted. 
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Appendix  B 
Tra.iectory  Matching 


Calculation  of  AV-^  and  AV2 

The  performance  capability  of  each  stage  is  calculated 
using  the  ideal  velocity  equation,  with  a correction  for  fi- 
nite bum  losses.  The  ideal  velocity  equation  ist 


AV  = Isp  ^0  ZTl 


(2-1) 


where  I is  the  specific  impulse*  g0  is  the  gravitational 
constant*  mQ  is  the  total  mass  before  ignition*  and  m^  is 
the  total  mass  after  burnout.  If  ST  represents  vehicle  struc- 
ture mass,  PROP  represents  propellant  mass,  and  PL  the  mass 
of  the  payload,  then  the  ideal  velocity  equation  for  each 
stage  becomes: 


*V1  = ^pl 


(ST1+PROP1 
ST^+ST 


+st2+prop2+pl 


2+PR0P2+PL 


) 


(B-l) 


4V2  * Isp2  «0  *n  ( 


st2+prop2+pl 

st2+pl 


) 


where  all  the  values  on  the  right  side  of  Equation  (B-l)  are 
specified . 

The  actual  velocity  change  capability  of  a stage  is  some- 
thing less  than  the  ideal  AV,  due  to  gravitational  effects  dur- 
ing the  finite  bum  time.  To  account  for  finite  bum  losses, 
trajectories  resulting  from  a AV  applied  impulsively,  were  com- 
pared to  actual  integrated  trajectories.  By  comparing  the 
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energies  of  each  trajectory,  an  estimate  was  made  of  the  AV 
loss  due  to  the  finite  bums.  Using  this  technique,  the  loss 
for  a AV1  applied  in  a 160  nm  parking  orbit  was  estimated  to 
be  about  .1#.  This  correction  is  then  applied  to  AV.^  before 
impulsive  targeting  is  accomplished.  The  losses  during  AV2 
sire  negligible  because  of  the  higher  altitude. 

It  should  be  emphasized  that  this  correction  for  finite 
loss  is  only  applied  to  help  the  impulsive  targeting  generate 
the  best  possible  guesses  for  the  values  of  the  mission  var- 
iables to  be  given  to  NS01A. 

For  specified  thrust  magnitudes,  the  bum  times  (de- 
noted by  t^a  and  t^)  of  each  stage  are  fixed  and  are  ob- 
tained through  the  mass  flow  rate  relationships  ast 


m,  = 


(B-2) 


^ba 

*bb 


m 


01 


- m 


fl 


mn 


m02  ~ mf2 

III  r\ 


(B-3) 


where  m^  and  m2  are  the  mass  flow  rates  of  each  engine » and 
c^  and  c2  are  the  effective  exhaust  velocities  as  given  by 


C1  = Ispl  «0 
c2  ’ Isp2  60 


(B-4) 
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Limits  of  First  Bum  Plane  Change 

Using  the  value  of  AV1  corrected  for  estimated  finite 
bum  losses,  the  maximum  attainable  first  bum  plane  change 
is  one  followed  by  a minimum  energy  (Hohmann)  transfer.  That 
is,  using  the  minimum  component  (of  the  total  AV^  available) 
for  completing  the  transfer,  leaves  the  maximum  component  of 
AV^  available  for  executing  the  plane  change.  If  rep- 

resents the  maximum  possible  plane  change  by  bum  one,  its 
relationship  to  AV^  is  shown  in  Figure  B-l,  where  the  sub- 
script H denotes  the  Hohmann  value. 

Given  r1G,  r2C,  V1C  and  V2C,  computation  of  V1TH  is  as 


follows: 


aTH  * 


rlT  * r2T 


Figure  B-l.  Maximum  First  Bum  Plane  Change 
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1TH 
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2aTH 
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IT 


TH 


(B-6) 


(B-7) 


where  aTH  and  Em  are  the  transfer  orbit  semi-major  axis  and 
energy,  respectively*!  Then  using  the  law  of  cosines: 


2m  ax 


= cos 


-1 


r 


1C 


+ v 


1TH 


2 V1C  V1TH 


liifl 
H J 


(B-8) 


In  calculating  lifgmax  n0  ef^or't  i®  made  in  trying  to  find 
an  actual  value  which  takes  into  account  the  magnitude  of  AV2 
and  trajectory  matching.  Thus,  il/gmax  strictly  311  upper 
limit,  and  the  range  of  possible  transfer  trajectories  can 
be  no  more  than  for  0°  < <J/2  < ^2max*  This  "then  is  the  range 
that  is  sampled  at  1°  intervals. 

As  a note  of  interest,  the  upper  limit  of  the  overall 
plane  change  possible  using  both  burns  may  be  determined 
from: 


r2TH  ° / r^% 


ax 


C08 


V ^ ♦ 
-1  V2C 

2 

m 


2 A y 2 . y 2 

V2TH_____2_ 


V2C  V2TH 


] 


(B-9) 


(B-10) 


where  '|f4max  is  the  maximum  plane  change  possible  by  bum  two 
(AV2) . Given  the  altitudes  of  the  parking  and  mission  orbiti 
and  the  available  AV1  and  dV2,  the  total  plane  change  possi- 
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Here  the  governing  equations  used  in  targeting  the  trans- 


fer for  the  impulsive  &V  approximation  are  developed.  Num- 
erous spherical  trigonometric  relationships  are  necessary, 
and  their  verification  is  left  to  the  reader  in  order  that 
the  logical  steps  remain  uncluttered  by  extraneous  explana- 
tion. Similarly,  the  reader  is  referred  to  the  listing  in 
the  front  of  the  book  and  the  various  Figures  for  an  expla- 
nation of  any  unfamiliar  notation. 

The  relationships  between  vectors  and  angles  for  the 
first  bum  are  as  shown  in  Figure  B-2,  all  referenced  to 


the  local  frame. 
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a 


b 


PT 


eT 


cos-1 

( C08  ij/g  COS 

(B-ll) 

sin-1 

[-Vic  “l] 

(B-12) 

rr  - a 

- P1 

(B-13) 

J^2 

+ V1C2  - 2V1C  &V1  cos  b 

(B-14) 

CM 

Eh 

iH  CM 
> 

_ _U_ 
rlT 

(B-15) 

2% 

(B-16) 

riT  v: 

1T  cos  lix-L 

(B-17) 

hT2 

(B-18) 

7i  - 

( Py  / Svp ) 

(B-19) 

Check  to  assure  a valid  transfer  orbit « 


1 


(B-20) 


If  both  conditions  are  satisfied,  continue.  Check  for  valid 
second  burn  conditions  (see  Pigure  B-6)i 


w = 


cos 


-1 


♦ V2C  - 

"5V 7 — 

ZV2T  V2C 


(B-21) 
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If  the  absolute  value  of  the  argument  of  Equation  (B-21) 
is  less  than  one,  continue. 

Notes  Other  quite  numerous  checks  must  be  made  through- 
out these  computations,  but  are  not  shown  here  (mainly  check- 
ing the  arguments  of  arc-cosine  and  arc-sine  terms,  and  quan- 
tities under  square  roots). 

Now  compute  the  true  anomalies  and  central  angle  of  the 
transfer  orbits 


/ / rii>)  ~ 1 

V1  = J 

j (pip  / r2<r)  ” ^ 
2 V e<p 


(B-22) 


°a  - v2  ' V1 


(B-23) 


Now  compute  the  corresponding  plane  change  required 
(tjj^)  at  the  second  burn  (see  Figure  B-3),  so  that  the  total 
two  burn  plane  change  is  equal  to  9 . 


--1  rrr»j 


(B-24) 


c = 


J* cos  [1/2  (Ca  - d)]  n 

2 tan  Lcos  Li'/T(Ca  Td)]  tan  [l/T  <9  + 


^4 


= TT  - C 


2. 

(B-25) 

(B-26) 


The  second  bum  flight  path  angle  is  obtained  from 


i|/-  = cos 


1 [rVl 

L2T  v2T  J 


(B-27) 
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Figure  B-3.  Transfer  Geometry 

Here,  the  angles  ^ and  the  vectors  V2T,  V2C  and  &V2  are 
all  referenced  to  the  e frame  as  shown  in  Figure  B-4.  The  e 
frame  is  a rotating  frame  that  describes  the  position  of  the 
IUS  at  the  second  burn  time.  The  er  axis  is  always  in  the 
r2^,  direction.  The  axis  is  in  the  plane  of  the  transfer 
trajectory,  perpendicular  to  r2T.  The  plane  change  angle 
(tjr^)  is  measured  in  the  ez  plane  from  l0.  The  second  bum 
flight  path  angle  is  measured  in  the  §r  $e  plane  from  eQ 

These  angles  are  shown  in  Figure  B-4. 

The  relationship  of  Equation  (3-1), 


AVgfAVi*  ijf^,  tj;2) 


(3-D 


is  contained  within  the  preceding  oribtal  transfer  equations. 
Its  solution  constitutes  the  (impulsive)  trajectory  matching 
process.  To  verify  this  relationship  and  gain  further  in- 
sight into  its  nature,  a closed  form  solution  for  the  copla- 
nar  case  is  derived  later  in  this  Appendix  (as  Equation  B-57). 

For  a specified  value  of  tJ/2*  the  value  of  tj/^  which  sat- 
isfies (3-1)  is  found  by  accomplishing  a one -dimensional  pa- 
rameter search  on  AVg  by  varying  . Rather  large  steps  are 
taken  until  the  solution  is  bracketed  by  two  values  of 
Then  the  final  solution  is  obtained  by  using  a Regula-Falsi 
(secant)  numerical  algorithm  (Ref  4i178). 

The  parameter  search  is  accomplished  by  using  the  re- 
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lationships  of  (B-ll)  through  (B-27)  to  determine  values  for 
the  angles  ijj^  and  t]/^  in  the  6 frame,  which  are  generated  by 
each  value  of  i|/^  during  the  one -dimensional  search. 

The  angles  tyo  and  establish  the  vectors  V2C  and  V2T 
in  the  § frame.  Thus,  the  value  of  AV2  corresponding  to  that 
value  of  ijj^  is  obtained  from  the  magnitude  of 

W-2=V2C  - V2T  (B-28) 

Once  the  one -dimensional  search  has  roughly  bracketed 
the  solution  to  Equation  (2-2),  then  the  Regula-Falsi  itera- 
tions are  used  until 

I *V2  k+1  " AW2*\  < ACCURACY  (B-29) 

where  AV2  is  the  actual  value  of  the  second  stage  velocity 
change  capability,  and  AV2  k+1  is  the  iterated  value  from 
the  Regula-Falsi  algorithm,  which  is 

Vl  - +*  - [fT2^|--  ] (4V2k  - 4V2‘> 

When  Equation  (B-29)  has  been  satisfied,  the  trajectory 
match  has  been  accomplished  and  the  values  of  the  angles  tj/^, 
t|/2»  and  are  all  known. 

The  last  step  is  to  solve  for  the  thrust  direction  angles. 
For  the  impulsive  AV  approximation,  the  thrust  directions  are 
the  same  as  the  velocity  change  directions  for  each  burn. 

The  components  of  AV1  (expressed  in  the  XYZ^  frame)  are 
obtained  from 
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*V1  - V1T  - V1C 


(B-31) 


where  V1C  = V1C  j j and  the  components  of  V1T  are  given  by 


V1TX  = V1T  sln  ^1 
V1TY  = V1T  008  ^1  cos  ^2 


^1TZ  = ^1T  sin  i \i f 


(B-32) 


The  direction  of  is  the  same  as  that  of  in  Fig- 
ure 4- 2,  so  that  the  first  bum  thrust  direction  angles  are 


found  from 


Cp2  = sin"1  (AV1Z  / AVX) 

= sin"1  (AV^y  / cos  $2) 


(B-33) 


In  analogous  fashion,  the  components  of  AV2  expressed 
in  the  6 frame  may  be  found  from  Equation  (B-28),  and  the 
known  values  of  tjj^  and  Then  a coordinate  transformation 

takes  w?  to  the  local  frame  through 


XYZ,  XYZ^  g 

2 = Ce  *V 2 


(B-34) 


The  coordinate  transformation  matrix  is  given  by 


XYZ  t - 


c y2  c YX 


1 3 8 Y2  c Yi 

U Yj. 


-s  y2  -c  Y2  s Yl 

c Y2  -s  Y2  8 YX 


(B-35) 


where  cosine  is  abbreviated  as  c,  and  sine  as  si  and  the  angles 


Y^  and  Yg  are  as  shown  in  Figure  B-4. 

Manipulation  is  necessary  to  obtain  the  angles  and 
Y2*  a top  and  side  view  of  Figure  B-4  is  shown  in  Figure  B-5. 
Here  a is  the  angle  in  the  plane  of  the  transfer  trajectory 
between  r2C  and  the  Y^Z^  plane,  so  that 


a = Ca  -.w  / 2 (B-36) 

r2TP  = r2T  cos  a (B-37) 

r gg  = r2rpp  sm  ij J2  (B- 38) 

Y = sin  ( r 2z  / r2ip)  (B—39) 

r 2X y = r2T  (B— 40) 

r2Y  = r2TP  C0S  ^2  (B-41) 

Y gp  = cos  ( r 2y  ^ r2XY^  (B— 42) 

Y2  = n / 2 + Y2p  (B-43) 


Once  AV2  is  expressed  in  the  local  frame,  the  second  bum 
thrust  direction  angles  are  obtained  fromt 

CP4  = sin'  (aV2Z  / ^2) 

(B-44) 

Cp-j  = sin-1  (aV2Y  / ^V2  008  $4) 


Checking  for  ambiguity,  if  cp-j  < 0,  then  cp^  = -(tt  + cp^) . 

Finally,  the  transfer  time  of  flight,  TOF,  is  found  from: 


'cl 


cos 


-1  |~eT  + 008  V1  ] 
Ll  + eT  cos  vx  J 


(B-45) 
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The  Functional  Relationship  AV^Ci^) 

The  purpose  of  this  section  is  to  yield  further  insight 
into  the  trajectory  matching  concept  by  developing  a closed 
form  relationship  for  AV2(^)  for  the  coplanar  case,  and  dem- 
onstrating that  a similar  relationship  holds  for  transfers  be- 
tween non-coplanar  orbits  when  the  amount  of  first  bum  plane 
change  (i|;2)  is  fixed. 

Figure  B-6  shows  the  three-dimensional  relationships  be- 
tween the  associated  vectors  and  angles  at  each  burn.  For  a 
coplanar  transfer,  angle  reduces  to  tj;^,  and  angle  w reduc- 
es to  The  known  quantities  are  r1T,  r2T,  V1C,  V2C,  AV1 

and  AVg. 

The  velocities  are  related  thru  the  law  of  cosines  (for 
coplanar  transfers)  as 


AV  ^ = V ^ + v 
avl  V1T  v] 


- 2 V 'T 
* IT  1C 


cos  tJi, 


(B-4?) 


AV  ^ = V ^ + V ' 
V2  V2T  V2C 


— 2 V2T  V 2C 


(B-48) 


Figure  B-6.  Three-Dimensional  Vector  Relationships 


The  total  energy  and  angular  momentum  of  the  transfer 
orbit  may  be  written  as 


V ‘ 
v IT 

hT  = 2 


IT 


2T 

2 


2T 


hT  = rlT  V1T  C0S  = r?T  V9T  008  ^ 


2T  2T 


(B-49) 


(B-50) 


The  task  is  to  develop  an  analytical  relationship  between 
AV2  and  ijjlf  with  as  the  only  unknown.  Rewriting  Equation 
(B-47)  as 


; 

m 

m- 


V,m2  - 2 V1T  Vir  cos  ij;-,  = AV.,2  - V1C2  (B-51) 


IT 


and  completing  the  squares  gives 


■ 

W 


(V1T  - V1C  cos  ^)Z  = AVX2  - 


V1C  + V1C*  cos*  (B-52) 


Now  taking  the  square  root  of  both  sides,  using  geometrical 
arguments  to  determine  the  proper  choice  of  signs,  gives  V1T 
as  a function  of  tjj^  onlyi 

V1T  = V1C  008  ^1  + ‘ V1C2  + V1C2  cos2  ^1  (B“53) 

o 

Using  the  energy  equation,  V2T  may  be  expressed  as 


V ^ - is.  U _ iLli  + y ‘ 
V2T  r2T  r1T  + V1T 


(B-54) 


and  solving  for  cos  tjj^  from  the  angular  momentum  equation 
give  8 


r1T  viT  c°s  ij;-, 


cos  *3  - V 


2T  2T 


(B-55) 
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Substitution  of  Equation  (B-55)  into  (B-48)  yields 


AV  2 = V 2 + v 2 
V2T  V2C 


2 V2C  rlT  V1T  C0S  ^3 

^ om 


(B-56) 


Finally,  Equations  (B-53)  and  (B-54)  are  used  in  (B-56) 
to  produce  (after  tedious  manipulation)  the  desired  result  of 
AV2  as  a function  of  ^ only« 

av2(^1)  = 2 cos  4i1  CV1C2  r2T  cos  \|»x 

~ V1C  V2C  rlT  + •Avi2  - V1C2  + V1C2  cos  11/^ 

^V1C  r2T  " V2C  riT^/r2T-^  " ^2  ^/rlT^ 


+ AV, 


(B-57) 


Equation  (B-57)  is  a closed  form  relationship  between 
the  first  burn  flight  path  angle  and  the  second  burn  velocity 
change  for  coplanar  orbits. 

An  analogous  procedure  can  be  used  to  derive  the  form 
of  AVgtij^)  for  non-coplanar  transfers.  The  major  difference 
is  that  now  Equations  (B-47)  and  (B-48)  become  the  three  di- 
mensional relationships 


*V12  = V1T2  + V1C2  " 2 V1T  V1C  008  *1 


(B-58) 


AV22  = V2T2  + V2C2  ' 2 V*r  V*r.  c°8  w ( B— 59) 


2T  2C 


where  B^  and  w are  expressed  by  (see  Figures  B-2  and  B-7) 
» cos-1  (cos  (Jfg  cos  ) (B' 
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Figure  B-7 . Second  Bum  Relationships 


These  relationships  show  that  now  i{i2  and  are  introduced  as 
additional  variables. 

Thus,  the  three-dimensional  equivalent  to  Equation  (B-57) 
is  a function  of  ij/2  and  in  addition  to  ij;^,  giving  it  the 
extra  degrees  of  freedom  that  allow  a range  of  solutions.  How 
ever,  careful  analysis  shows  that  if  the  value  of  ijj2  is  speci- 
fied, then  the  value  of  tj is  fixed  through  Equations  (B-24, 
B-25  and  B-26).  So,  by  specifying  the  amount  of  plane  change 
to  be  accomplished  during  the  first  bum,  the  amount  left  for 
the  second  bum  is  fixed  by  the  overall  plane  change  require- 
ment! and  again  AV2  becomes  a function  of  iju  only. 


•”  ■*  ■“ 


r 


Appendix  C 

Finite  Bum  Dynamics 


Equations  of  Motion 

In  this  section  the  nonlinear  vector  differential  equation 
of  motion  is  developed  and  then  put  into  a state  variable  for- 
mat for  implementation  of  the  numerical  integration  routine 
(SET/STEP).  The  local  (XYZ^)  frame  is  used  as  the  inertial 
reference.  In  this  frame,  the  position  vector  and  thrust  vec- 
tor are 


r(t)  = rx  i + ry  j + r2  k 
T = Tx  I * Ty  3 + Tz  k 

Thus,  finite  bum  dynamics  are  described  by 
Fg[r(t)]  + T = m(t)  r(t) 


(C-l) 

(C-2) 


(C-3) 


where  m(t)  represents  the  mass  of  the  vehicle  at  any  time,  tj 
and  Fg  is  the  force  due  to  gravity,  expressible  as 

F = TiLJlLj?  r(t)  (C-4) 

g r(t)-* 

where  ix  is  the  gravitational  parameter. 

Since  T is  constant,  the  mass  flow  rate  is  constant, 
yielding 

m(t)  * mQ  - mt 

Thus,  Equation  (C-3)  may  be  expressed  as  (omitting  time  de- 
pendency) 
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r = ziLJ?  + T 

r.3  nu  - rot 


(C-5) 


Expressing  (C-5)  in  component  form  gives 


r3  rx  + mQ  - At 


^ ry  + 


(C-6) 


~3  rz  + itu  - At 


Finally,  using  a six  element  state  vector  (X)  where  the 
first  three  elements  are  the  position  components,  and  the  last 
three  elements  are  the  velocity  components,  the  equations  of 
motion  become t 


X1  " X4 
X2  = X5 
x3  = x6 


-4  X, 


4 " r3  mQ  - mt 

-M,  T 

*5  r3  m0  - At 


-4  X. 


(C-7) 


,3  mn  - At 


During  the  coast  phase  of  the  trajectory,  the  components 
of  T are  all  zero. 

The  components  of  the  thrust  vector  are  calculated  from 
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the  current  values  of  the  thrust  direction  angles  as  follows: 
For  the  first  bum 


Tx  = T^  cos  rp2  cos 

Ty  = T1  cos  cp2  sin  (C-8) 

Tz  = Tx  sin  cp2 

and  during  the  second  bum  as 

Tx  = T2  cos  cp^  cos  ^ 

Ty  = T2  cos  cp^  sin  cp^  (C-9) 

Tz  = T2  sin  <*>4 

Thus,  the  means  of  accomplishing  the  targeting  for  the 
finite  bum  case  are  contained  in  the  equations  of  motion 
(C-7).  By  adjusting  the  thrust  directions  through  Equations 
(C-8)  and  (C-9) • and  the  transfer  TOF  by  varying  the  integra- 
tion time,  NS01A  repeatedly  demands  integration  of  Equations 
(C-7)  until  values  for  the  mission  variables  are  found  that 
produce  the  exact  mission  orbit  desired. 

Transigr  Angie 

The  transfer  angle  (TA)  is  defined  as  the  angle  between 
the  vectors  ric(t)  and  r2C(t  ♦ TOF),  where  finite  bum  dynam- 
ics apply.  Thus,  TA  is  slightly  different  from  the  central 
angle  (C  ),  used  to  describe  the  angle  between  these  two  vec- 

a 

tors  during  impulisve  targeting. 

To  find  the  value  of  the  transfer  angle  resulting  from 
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targeting  step  one  (where  it  was  left  free),  a third  vector 
(Sr)  is  first  formed,  as  shown  in  Figure  C-l.  In  the  local 
frame,  the  component  form  of  r-^t^)  is  always  just  rlc  T. 

The  components  of  + TOF)  are  known  after  numerical 

integration  of  the  transfer  trajectory. 

The  vector  Ar  is  then  determined  from 

^ = ?2C  " ?1C  (C-10) 


and  the  transfer  angle  is  obtained  using  the  law  of  cosines  as 


cos 


-1 


r 2 2 

r2C  rlC  ~ 
2 r1G  r2C 


(C-ll) 


During  the  second  step  of  targeting,  if  the  transfer  is 
between  coplanar  orbits,  it  is  necessary  to  compute  the  trans- 
fer angle  as  a function  of  time,  expressed  in  the  perifocal 
(PQW)  frame.  In  this  case.  Equation  (C-ll)  takes  the  form 
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TA(t)  = cos 


(C-12) 


-1 


r2C2(t  + TOF*)  + rlc2(t)  - Ar2(t) 
2 rlc(t)  r2C(t  + TOF*) 


and  is  solved  for  the  first  bum  time  (tfcl)  such  that 

TA(tbl)  = TA*  (C-13) 

Computation  of  the  vectors  ric(t)  and  r2C(t  + TOF  ) is  de- 
scribed in  Appendix  D. 


Appendix  D 
Targeting  Step  Two 

Phasing 

Phasing  is  simply  a matter  of  selecting  the  proper  launch 
position  ric(ttl),  such  that  the  trajectory  targeted  in  step 
one,  when  initiated  at  that  position,  satisfies  either  the  con- 
straint n2,  or  the  constraint  TA(t).  Since  r-^t)  is  known, 
this  becomes  a problem  of  finding  the  proper  time  (ttl)  for 
first  stage  ignition. 

The  sequence  of  operations  carried  out  in  phasing  (and 
their  significance)  are  listed  here  in  the  order  that  they 
are  accomplished t 

1.  From  the  given  orbital  elements  of  each  orbit,  ( "t ) 
and  r2C(t)  (when  needed)  are  computed  in  their  respective 
perifocal  (PQW)  frames.  The  vector  ric(t)  tracks  the  IUS 
in  the  parking  orbit,  and  ?2C(t)  tracks  the  mission  orbit 
target  position,  if  specified. 

2.  A value  for  t^^  is  found  such  that  becomes 

the  proper  launch  point  to  satisfy  either  n2  (if  the  or- 
bits are  non-coplanar) , or  TA*(t^)  if  transferring  be- 
tween coplanar  orbits.  In  addition  to  specifying  the 
mission  variable  t^,  this  step  fixes  the  orientation 
between  the  local  frame  and  the  geocentric-equatorial 
frame  (see  Figure  D-l) . 

3.  Through  knowledge  of  the  orientation  between  the 
frames,  a coordinate  transformation  is  accomplished  to 
convert  the  thrust  direction  angles  expressed  in  the 
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Figure  D-l.  Orientation  Between  Frames 

local  frame,  to  their  corresponding  values  in  the  geo- 
centric-equatorial frame. 

The  actual  accomplishment  of  the  phasing  steps  listed 
above  involves  a great  deal  of  computations,  and  only  the  high- 
lights wifl.1  be  given  here.  The  reader  is  referred  to  the  com- 
puter code  listing  if  more  details  are  desired. 

To  compute  ric ( t ) and  r2C(t)  in  their  perifocal  (PQW) 
frames,  epoch  time  (tQ)  is  used  to  fix  their  initial  positions. 
At  epoch,  rlc(tQ)  is  directed  along  the  P^  axis,  and  r2C(tQ)  is 
at  an  angle  u02  with  the  ?2  axis.  First,  the  angular  veloci- 
ties are  determined  from 
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t 


(D-l) 


(!)-,  = 


lie 

•ic 


(l)o  = 


2C 

2C 


Then  the  two  vectors  can  be  calculated  at  any  time  t as 
rlC(t)  = (rlC  cos  (1)^  t)  sin  t)  (D-2) 

r2C(t)  = [r2C  cos  (d)2  tf  + uQ2)]  P2 

(D-3) 

+ [r2C  sin  (w2  tf  + u02)]  Q2 

* # 

where  t-  * t + TOF  » and  TOF  is  the  transfer  time  of  flight 
as  targeted  in  step  one. 

If  n2  is  to  be  satisfied,  then  tfcl  may  be  found  as  fol- 
lows i 


(D-4) 
(D-5) 

Here  u1  represents  the  angle  that  ?ic(t)  generates  in  going 
from  tQ  to  t^i  and  is  the  value  resulting  from  step  one, 
where  longitude  of  the  ascending  node  was  left  free. 

If  TA(t)  is  to  be  satisfied  (coplanar  orbits).  Equation 
(C-12)  is  used.  A one -dimensional  search  on  t is  made  until 
a solution  is  bracketed*  then  a secant  technique  is  used  to 
get  the  exact  t^  such  that  Equation  (C-13)  is  satisfied. 

The  last  step  is  to  transform  the  thrust  angles  from 
the  local  frame  to  the  geocentric-equatorial  frame.  Once  ttl 


ul  = n2  - 
% 

*bl  “ 
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is  known,  the  orientation  between  the  frames  is  fixed  by  the 


; 


given  orbital  elements  of  each  frame,  and  ) as  shown 

in  Figure  D-l.  Here  v^  is  the  single  between  the  P1  axis  and 

rlC  ^ ^bl ^ ’ 

To  make  the  transformation,  unit  thrust  vectors,  uT1  and 
uT2,  are  formed  (expressed  in  the  local  frame)  from  the  tar- 
geted thrust  angles.  Then  these  vectors  are  each  transformed 
to  the  geocentric-equatorial  frame  by 

_XYZ  XYZ  _XYZ£ 

UT  = CXYZX  UT  ^D_6) 

After  the  transformation,  angles  are  formed  again  from  the  unit 
vectors. 

The  three  Euler  angles  between  the  frames  are  (by  order 
of  removal),  v^,  i^  and  The  angle  v-^  is  obtained,  once 

t^  is  known,  from  the  relationship 

V1  = wl  Hi  (D-7) 

The  transformation  matrix  is 

(cfl1cv1+s01ci1sv1)  (cn1sv1-sn1ci1cv1)  (3^3^) 
(sn1cv^-c01ci1sv1)  (sC^sv-j+cf^ci-LCv-^)  (-si1cfi1) 
(-si1sv^)  (sijcv^  (ci-^ 

(D-8) 


XYZ 

CXYZt  = 


Mission  Delay  Retargeting 

For  transfer  between  non-coplanar  orbits,  the  mission  op- 
portunity times  occur  at  equal  intervals  after  the  first  ttl 
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time.  This  is  because  the  first  bum  must  occur  at  the  same 
position  in  the  parking  orbit  each  time.  The  time  interval 
between  the  possible  launch  times  is  just  equal  to  the  period 
of  the  parking  orbit,  which  is  given  by 


(D-9) 


and  the  thrust  direct  angles  are  the  same  at  each  possible 
launch  time. 

For  transfer  between  coplanar  orbits  (rendezvous  assumed), 
the  mission  opportunity  times  are  again  periodic,  but  now  the 
interval  is  equal  to  the  synodic  period,  which  is  found  by 


p — 2 TT 

s - |0>1  - m2| 


(D-10) 


In  this  case  the  launch  position  changes  each  time,  so  that 
for  each  sequential  opportunity,  the  corresponding  angle  v^ 
must  be  computed,  and  the  coordinate  transformation  of  the 
thrust  angles  repeated. 

Summarizing,  the  mission  opportunity  times  (times  when 
proper  phasing  occurs)  are,  for  non-coplanar  orbits 


= K * p 


(D-ll) 


and  for  coplanar  orbits 


tv,  = tv,  + PH 
blk+l  blk  8 


(D-12) 


fox*  lc  = If  2,  • • • • 
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Appendix  E 

Accuracy  Determination 

Thrust  Misalignment  Sensitivities 

The  nominal  values  of  the  thrust  direction  angles  (p^, 

P2.  P^  and  $4)*  expressed  in  the  local  frame,  are  obtained 
in  targeting  step  one.  The  procedure  used  to  obtain  inser- 
tion error  sensitivities,  due  to  thrust  vector  misalignment, 
is  to  perturb  these  nominal  values  slightly  by  assuming  a one 
milliradian  (mr)  alignment  error  about  each  axis  of  the  XYZ^ 
f rame , in  turn . 

To  generate  perturbed  thrust  angles,  nominal  unit  thrust 
vectors  (u,^  and  uT2)  are  first  formed,  with  their  components 
expressed  in  the  XYZ?  frame.  Then  these  vectors  are  coordi- 
natized  in  a frame  which  is  misaligned  by  one  mr  about  the  X£ 
axis.  Next,  their  components  in  the  misaligned  frame  are  con- 
verted to  misaligned  thrust  angles.  Lastly,  Equations  ( C -7 ) 
are  integrated  using  these  misaligned  thrust  angles,  so  that 
the  insertion  error  sensitivities  due  to  a one  mr  misalign- 
ment about  the  axis  are  obtained.  This  gives  the  first 
column  of  the  sensitivity  matrix  in  Equation  (5-2).  The  same 
process  is  then  repeated  for  misalignments  about  the  Y^  axis 
and  the  Z ^ axis,  generating  the  second  and  third  columns  of 
the  matrix. 

The  insertion  error  sensitivity  matrix  of  Equation  (5-2) 
is  then  partitioned  as  shown,  to  form  the  individual  position 
and  velocity  matrices.  Next,  the  maximum  eigenvalues  of  these 
matrices  are  determined  by  use  of  a standard  subroutine  designed 
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for  that  purpose,  and  the  two  "worst  case"  sensitivities  are 
formed  using  Equation  (5-7) • The  derivation  of  this  equation 
will  now  be  explained. 


The  derivation  is  motivated  by  formulating  the  problem 
as  Equation  (5-6) 


y = A x 


(5-6) 


where,  for  a fixed  magnitude  of  the  vector  x (thrust  misalign- 
ment vector),  it  is  desired  to  find  the  maximum  possible  mag- 
nitude of  the  vector  y (representing  insertion  position  or 
velocity  error),  caused  by  the  sensitivity  matrix,  A. 

To  express  the  worst  case  sensitivities  with  respect  to 
milliradians  of  misalignment,  the  problem  becomesi  Given 


— / 2 2 2 

lx  i = Vx1  + x2  + x3  = 1 mr 


(E-l) 


find  the  associated  value  of  iyimax* 

The  steps  in  the  development  of  Equation  (5-7)  are  as 
follows 


y = A x 


(5-6) 


_rp_  fp  rp  ... 

y y = x A A x 


(E-2) 


thus  AXA  is  a symmetric  matrix  where 


ata  = m_1  a m 


(E-3) 


where  A is  the  diagonal  matrix  whose  elements  are  the  eigen- 
values of  ATA,  and  M is  the  normalized  (and  dimensionless) 


■ 


m 


modal  matrix  of  AaAj  and  is  scaled  such  that 


M"1  = MT 


(E-4) 


I 


substituting  (E-3)  into  (E-2)  gives 


= ? II'1  A « ; 


now  let 


z = M x 


(E-5) 


(E-6) 


f, 

i' 


I 


and  taking  the  transpose  of  each  side 


- i*  MT  = J1  M-1 


(E-7) 


Now,  substituting  z and  z1  into  (E-5)  gives 


— T—  — 2 — T — 

y y = i y i = z A z 


(E-8) 


and  performing  the  vector  matrix  multiplications  yields 

y2  = ^ z2  + X2  z22  + Zz  (E-9) 


I 


But  the  vectors  x and  z are  related  through 


z^z  = ar1'  (NITM)  x = ^x 


since  M is  an  orthogonal  matrix  such  that 


MTM  = M-1  M = I 


where  I is  the  identity  matrix.  Thus 


l z l * ixi 


(E-10) 


(E-ll) 


(E-12) 
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The  solution  involves  checking  the  four  possible  cases,  which 
are  i 

Case  1«  z-^  = 0,  z2  = 0,  z^  = 0.  This  case  is  impossible 

since  it  violates  the  constraint  Equation  (E-14). 


Thus,  the  maximum  magnitude  of  y must  be  determined  by  the  larg- 

T 

est  of  the  three  eigenvalues  of  the  A A matrix,  and  is  always 


iyi  = vT 


(E-23) 


for  an  Ixi  of  one  mr. 


Thrust  Magnitude  Sensitivities 

Here  the  magnitudes  of  the  thrust  vectors  of  each  stage 
( and  Tg)  are  changed  by  ±1#,  and  Equations  (C-7)  integrated 
for  the  nominal  time  (TOF  ) so  that  insertion  errors  may  again 
be  computed.  The  thrusts  of  both  engines  are  assumed  to  devi- 
ate by  the  same  plus  amount  or  the  same  minus  amount,  so  that 
insertion  position  and  velocity  error  sensitivities  may  be  de- 
termined for  both  positive  and  negative  thrust  magnitude  devi- 
ations. 

The  bum  times  in  Equations  (C-7)  are  also  changed  by 
thrust  magnitude  variations  and  must  be  recomputed  prior  to 
each  integration  as 


A - II* 

lp  c. 


m - Zfc 

2p  ~ c. 


(E-24) 


m01  ~ mf] 


m02  " mi 


(E-25) 


where  the  subscript  p denotes  the  perturbed  value, 
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TDeH*PI*S3Rt  <ATH**3/ENU> 

COMPUTE  0*1  1N3  0*?  (IOEAD,  AND  THE  SRM  9 JRN  TIMES* 

»R3»l«PROPt/30  S PR OP 2« PRO*? /GO  i PL«P./68 
Mt*5Tl*PR0»t 
M2*ST2*PR0»2 
M01«M1*M2*»L 
M02«H24PL 
MPl*M01-PR3Pl 
MF?*M02-PR0P2 
C1«ISP1*GP  • 

C2«13P2*G0 
M91«T1/Cl 
M3?*T2/C2 
T9A«(M01-NFll/MOl 
TB3» ( i02«Ms2l  /MO  2 
9V1*ISP1*39*1.3G (M01/HF1I 
9V?«ISP2*59*1L3G<M02/HF2> 

CORRECT  FOR  FINITE  BURN  LOSSES  OF  0*1*  FINITE  LOSSES  OF  0V2  NEGL. 

FILOSS*.O01 
9*1 *9*1 • (9*1* rI LOSS) 

COMPUTE  TNE  MAXIMUM  P0SSI9.E  PLANE  CRINGE  JURING  THE  FIRST  BURN* 

SY2MAX*AC3S((*lC**2*VlTH**?-9Vl**2)/(2.*Vi:*VlTH>) 

NY*  MAX*AC9S ((**C**2**2TH**?»3*2**2) /(2.**2C**2TH)) 
SYrot*SY2NlX*S*lNAX 

°RTNT  INPUTS! 

°RINT*,-XXXXXXFXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
3XXXXXXXXXXXXXXXXXXXXXXXXXXXXX1* 

°RINT*,“XXXXXXXXXXXXXXXXXXXXXXXXXXXXX  OUTPJT  XXXXXX<XXX<<XXXXXXXXX 
NXXXXXXXXXXXX ‘XXXXXXXXXXXXXXXX" 

PRINT* 

PRINT* 

:p(«»AXPL.E3.  0)30  TO  JN1 

®RTNT*,"TMIS  IN  A MAXIMUM  PAf.OAO  RUN  U5IN3  MIT  CONOITION  SET  TNO" 
PRINT* 

PRINT* 

041  CONTINUE 

PRINT*, "TME  3R9IT  PARAMETERS  ARE*  “ 

PRINT* 

PRINT*, “HI*  NM.  AN3  M?«  “,M2,-  NM. " 

®».INT**-Ri:*  -,*1C/CNMF,*  N,M.  ANO  R?C«  **,  R2C/CNNF, “ N* M«  * 

PRINT*, “*13«  ",V1C,"  FT/SEC.  ANO  *23*  -,*23,"  *T/SEC.  * 

PRf  NT*  ,**CO,LlN*  • MOMMANN  TRANSFER  V£L9*TTI-S  IRE!  " 
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NM* 


( 


o 


( 


! 1 

• I 

: i 


IFUUNP.GE.  t)G3  TO  123 
C 

C CALL  SU8R3 JTIN:  AGUESS  NHISH  C0«»UT£S  INITIAL  GUESS'S  FDR  NS01A 
C 

CA.L  AGUESS  C*HT1,PMI2, PHI  3,  P-»IL,TOF> 

C 

C CNE3<  IF  TUtSfEA  POSSIBLE  FOR  THIS  VALUE  OF  5Y2» 

C 

IF(3A.EQ.0.>S3  TO  1 

C 

TO*«TOFmi 

9F»1. 

ST«1. 0E5 
Y(l»«PHIl/ST 
V<2»xPHI2VSF 
Y(3)»PHI3VSF 
V<4|«PMIkVSF 
y<i»*TOF/sr 
123  CONTINUE 
C 

C CALL  THE  N3N-. INEAR  EQUATION  SOLVING  ROUTINE,  NSOIA 

C 

CA.L  N«5OlA<5,V,F,AJINV,l.E-G,l.Er,l.S-S,250,0,H> 

C 

ACNE2K*  < A85I R20 A-R2CPI ) /3NNF 
I*C  ACM£CK..T.1<’,>G0  TO  19 

•RINT*, “FOR  sr?«  ",SY2*R0,"  NO  FINITE  M»j?;T0RY  is  possible- 
print*, -ache:<*  acheck,*  n.h.  difference  * 

IFU'JHP.GT.OtSO  TO  30 
PRINT* 

PRINT* 

GO  TO  t 

19  CONTINUE 

To*«rc5>*sr 

T2J.TOF-TM 

®Hll«V<l»MF 

•HI2«V (2)*SF 

•HH«VI3>*SF 

PNIA«VIk)*SF 

•91  NT* 

•RINT* 

C5Y2»SY2*RJ 
PRINT ( 6,?kt  :ST2 


2k  roRNATIPX, ••••'••••  FOR  FIRST  BURN  PLANE  CNAN3E  OP»  -,F5.1,-  OEG 

i.  MtMHIMMM-,///) 


PRI N**  .“THE  PINAL  ACOUIREO  ORBIT  DATA  ISl  " 

•RINT* 

PRINT*, -FINA.  POSITION  AN J VELOCITY  ARE*  - 

PRINT*  , "POSITION  IN  XYZI  X (1WCN1P,X(  ?) /CNMF,  XC  3)  /5NNF 

PRINT*, -VELOII'Y  IN  X YZl  (M ,X<5» ,X(Sl 

PRINT*,**  R?C:«  -,R2CC/CNNp,-  N.N.  - 
PRINT*,-  V?C:»  -,V2CC»"  FT/SSC.  " 

•RINT*, -»?;a«  “ ,R2CA/CNNr • - N.N.  - 

•RIN'*,-P*;P*  -,R2CP/CMNF,“  N.N.  - 
•RINT*, -Il*,«  ",  II  A*RO,  “ )E  ••  * 

PRINT*, "I*;«  ",  1 2*R0, “ 3EJ.  " 

PRINT*, -023«  ", 02G*R0,-  3E3.  " 

PRINT*, “F2CC*  -,E2CC 
•RINT*, -TA«  ", TA*RO»“  0E5.  * 

PRINT* 


A. 


132 


I 


! i 


t 


•RtNT* 

PRINT* 

PRINT*, "THE  RESULTS  OF  TAR3ETING  ARE*  " 

PRINT* 

PRINT* 

°HI1 F*PHll*RO 
“■II2-«PHI2*R0 
®HI3F*PHI3*R3 
®HI AF«PHIA*R0 

PRINT*, "THE  FIRST  STAGE  THRUST  OIRECTION  ANGLES  ARE*  - 
PRINT* 

PRINT  (5,  751  »HI1F,»HI2F 

25  rORNAT  I10X,"*HI1*  ",F10.5,"  OEG>  ",S<,"»NI2*  "PlD.S,"  OEG.  “» 

®RI  'IT* 

PRINT* 

PRI NT*  t "THE  SECONO  STAGE  TNRJST  DIRECTION  AN3.ES  ARE*  “ 

PRUT* 

PRINT (6,65)  »HI3F,PMI5F 

65  e0RiATC10X»"*Hr3*  ",F10.J,"  DE3, ", SX,"PNI6*  "*16. S,"  OEG.  ") 

PRINT* 

PRINT* 

PRINT*, "BURN  MO  OCCURS  AT  TINE*  “,T22,"  SEC.  " 

PRINT* 

PRINT*, “TOTA.  r OF  IS*  ~»TOe,"  SEC.  “ 

PRINT* 

PRINT* 

C 

C CALL  SUBROUTINE  ERROR  WHICH  COMPUTES  INSERTION  ERROR  SENSITIVITIES 
C 

CALL  ERRORS ( ®HI 1 , PHIZ, PHI 3, PHI  A, TOP) 

PRINT* 

PRINT*, “THE  RESULTS  OF  THE  ERROR  STJOf  IS  THE  SENSITIVITY  MATRIX*" 
PRINT*, -UNITS  * RE  N.N.  OR  fT/SEC.  PER  MtLIRAOIAN  OF  MISALIGNMENT." 
PRINT* 

PRINT* 

DO  20  K*1,S 

PRINT (6,111 ( SNS ( K,  J ) , J*  1,  3) 

11  P3RM«T (10«,313<,F13.A),/) 

?0  CONTINUE 

PRINT* 

PRINT* 

PRINT*,"— ———————————————— 

C 

C NOV  COMPUTE  TM-  WORST  CASE  INSERTION  ERRORS  USING  EtGINVALUES 

C 0?  THE  SENSITIVITY  MATRIX 

C 

PRINT* 

PRINT* 

00  3 1*1,3 
00  k J*l,3 
A(T, Jt*SNS<I,J) 

<*I  *3 

0( I , J) *SNS  <X, J* 
h CONTINUE 

3 CONTINUE 

CALL  NTSPO  , 3,3,  ATP* 

CALL  MTSP(1,S,T,8TP) 

CALL  MMPYCAT*,A,C,3,3,3) 

CALL  NHPYI 3T*»  3, 0,3, 3, 3) 

▼NO»l 
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CA.l  RGEIS<3,I,:,INO,EIGVA.#EIGVEC,N3RO 
PRINT*, “FOR  3ELTA  R HISS!  " 

00  5 1*1,3 
PRINT* 

PRINT*, “EISINVAIUE  HO.  “,I»"  AND  ITS  ASSOCIATED  EIGXNVECTOR  IS«" 
PRINT* 

PRINT*, EIGVAim 
PRINT* 

PRINT*, (EIGVEC(J,I)  ,J«1,3> 

CONTINUE 

PRINT* 

AE1*CABS(EXGVAL(1)>  . 

AE?*5A*S(EIGV1.<2>> 

A£3*CA0S(£IGVA.(3n 

I- ( (AE1.GE. A?2» . AND.  I AEl.Gi. A£3> > EIGNAX-AE1 
Ir((AE?.GE.Am  .ANO.  (AE2.  GE.AE3))  EIGNAX*AE2 
T-( (AE3.GE. A El) . AND. ( AES. 3£. A £2) > EI3iA<*AE3 
0R1AX*SQRTCEISNAX> 

*RI NT*, "♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦ ♦♦*♦♦♦♦♦♦♦♦♦♦♦*»• 

PRINT* 

IN3*1 

CA.l  R3EI3(3,3» D,INO,EIGVA. ,£ISVEC,M3RK) 

PRINT*, "FOR  OElT  A V HISS»  " 

00  5 1*1.3 
PRINT* 

PRINT*, “EI5INVAIUE  NO.  “,I,"  AND  ITS  ASSOCIATED  EIGINVECTOR  ISA- 
PRINT* 

®RINT*,EI3VAl(TI 

PRINT* 

PRINT*, tEI3VE0IJ,I>  , J*l,3) 

CONTINUE 

“RINT* 

AEt  »3A0S (EIGVAl ( 1) > 

A£?«0A9S<EtGVAl<2>> 

AEHOAflSIEIGVA. C3») 

IF(<AE1.GE.AE»* . AND.(AE1.3E.1£3))EIGHAX*AE1 
T"<nE?.GE.AEl*  . AN3.UE2.G£.AE3)>SISiAX*AE2 
IF  ( ( AE3.  GE,  A III  • ANO.  (AE3.GE.l£2nEISNA(*A£3 
DViAX*SQ»T<£IGMAX> 

PRINT* 

®RI NT* 

■*RI  NT*  . — 

PRINT* 

PRINT* 

PRINT*, -Hi  43»ST  CASE  INSERTION  POSITION  ANO  VELOCITY  ERROR  - 
PRINT*, -SENSITIVITIES  OUE  TO  -THRUST  XIS#lIGN1£NT-  ERROR  AREA  “ 
PRINT* 

PRINT* 

PRINT* 

PRINT*, -ORNAX*  -,ORNAX,-  (N.H.)/(N,R.)  * 

PRINT*, "0V1AX*  ”,  3 VN AX  »" ' ( ?T/S£S. > f t N. R. » * 

PRINT* 

PRINT* 

PRINT*,” 

COMPUTE  TNRUST  MAS,  DEVIATION  ERROR  SENSITIVITIES 
VIA  SUBROUTINE  THRUST  FI 

CA.L  THRUSIF C R1C , VIC, TOF, Cl ,02) 
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C CONVERT  TARGETING  INFORMATION  TO  THE  3-5  PRANE 
c via  subroutine  burni 
c 

PRINT* 

PRINT* 

•RINT*,-TME  TARGETING  INFORMATION  CONVERTEO  TO  THE  GEOCENTRIC  " 
PRINT*  ,**£9JAT3RIAL  FRAME  f0R  THE  FIRST  FIVE  TRANSFER  START  TIMES*--- 
PRINT* 

®RIMT* 

IFCTMETA.3T. J.>30  TO  250 

CALL  BURNK  PHIl ( PHI2*  PHI3,PHIL,TOF,TA,  R1C,  VIC, R2C, V2C) 

GO  TO  251 

250  CALL  9URNN(PHI1,  PHI2, PHI3, »HI A.02G, RlC , V1C,R20,  V2CI 

251  CONTINUE 
PRINT* 

PRINT* 

•PINT* ,-*••••**••••***••*•*•••**••*•**•••**••*•••***•*•*•**•*•••** 

PRINT*, -*•••••••**••*♦♦•*•**•• ••••••••••*••**•*••••*•*•*••**•••••• 

J’JNP*  JUMP*t 
C 

C STEP  TO  THE  NEXT  INCREMENT  0*  FIRST  3 JRN  P.ANE  CMAM5E 

C 

GO  TO  1 

30  CONTINUE 
PRINT* 

PRINT* 

PRINT*, "THE  cJi.L  RANGE  OF  ’OSSIBLE  TRANSFERS  NAS  BEEN  EVALUATED,  - 
PRINT* 

PRINT* 

PRINT* ,**XXXXXXXXXXXXXXXXXX  ENO  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
5XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  " 

END 

SU3R0UTINE  A3UrSSIPHIl,PHI2,PNI3,PHH,,  TOPI 
C 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX<XXXXXX<XXXXX<XXXXXXXXXXXXXX 

C this  ROUTIN'  TARGETS  VALUES  r0R  THE  EMULSIVE  TRANSSER,  WHERE 
C FINITE  BURN  LOSSES  FOR  FIRST  STAGE  R JRN  ARE  INOLUOED, 

C XXXXXXXXXXXXXXXX XXX XX XXXXX<XXXXXXXXXXXXXXX<XXXXXXXX<XXX<<XXXXXXXXX 

c 

DIMENSION  <<S) 

CONNON/TT/EHJ,PI,CNMF,RO,OR,RlT,VlT,Ri:,Vi:,R»T,V2T,R2C,V2C,OVl, 
$0V2,TNUl,fNJ2,“»CA,AT,SVl,SY?,SY3,SV%, I2M, THETA 
REAL  I2N 
PRINT* 

PRINT* 

C 

C NON  CALL  SJ5R0JTINE  IGUESSl 

C 

5V1«0. 

SY3*0. 

CA.L  IGUESS 

C 

C CHECH  IF  TRANSFER  POSSIBLE  FOR  THIS  VALUE  OF  SY2I 

C 

IFCCA.EO.O.) 30  TO  269 

C 

C NON  COMPUTE  TOP* 

C 
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1 * 

|!  ( 

I j 


IF(E.  GT.  1.  >G3  to  4 

E:i*ACOSC<E*:3STTNUl>  >/<l.»£*COS(TNJim 
EC?*ACOS((EOOSCTNU2)  ) / Cl  • ♦E’COSCTNUZ ) ) ) 

T3r*SQRT(Ar  **T/ENU>* t Cfc32-E*SINCEC2>)  -{C31-E*SINCS;i> >> 

GO  TO  7 
C 

C FOR  HYPERBOLIC  ORBIT  TOFI 

c 

4 Yi*<E*cos<nm>  >/{i.+E*cos(T'Uin 

cl*  ALOGC Y1»S3RT CV1**2-1.)> 

Y?*CE*COSCTNJ2>  > /(l.  + E*C3S(TYJ2>» 
r?«AiOGCY2>S3RT (Y2**2-l.)> 

TOF*SQRT((-4n**3/EMU)*<CE*SIRH(e2)-F2)-li*SIYHCFl)-Fl»» 

7 CONTINUE 

C 

C THE  FIRST  3URN  THRUST  DIRECTION  ANGLES  ARE* 

C 

VI T X* VI T*SI N( SY 1 ) 

Vlf Y=V1T*C3S(SY1)*C0S  (SY2) 

V1TZ»V1T*03S(SY1)*SINCSY2) 

0V1X*V1TX 

•3VlV*ViTY-413 

0Y1Z*V1TZ 

0V1XY7*SQRT ( 0V1 X**2>0VlY**2OVlZ**2) 

*HI  2*  A SIN<  3 VI 7/0  VI) 

PHI1=ASIN(3V1Y/(0V1*C0S CFHI2) ) ) 

C 

C N34  COMPUTE  3V2  IN  THE  E FRAME  * 

C 

55  V2CX«0. 

V?CV*V2C*03SCSY4> 

V?CZ*V2C*SINCSV4) 

V2TX*V?T*StNCSY3) 

V2TY*V2T*C3S<SY3> 

v?rz=o. 

nv?XE»V2CX*V2TX 
0V2 YE=  V2CY*  V2TY 
0V?ZE«V2CZ-V>T7 

C 

C *134  OOOROHATI7E  0V2  IN  THE  INERTIAL  PRAM£I 

C 

-*i.«OA-PI/?. 

R2TS>=R2T*33S(3'L) 

”’?Z*R2TP*SIN(SY2) 

GAM A is  AS  INC R2Z/92T) 

’’YiRZTP^COSCSYE) 

R?XY*R2T»:3S<G4NA1> 

»?3*ACOS IR2Y/R2XY) 

■»?X«-R2XY»SIM(G?P) 

OV2XsOV2XE*C3S(GAHA2)*COS(3AMAl)-OV2Y‘*STN(GA4A2)-OV2?E*33SCGAMAZ) 
a**I NCGAHA1I 

TV’ Y*  B V2XE*SI N( G AHA2) •COSCCAMAl) ♦0V2YE*C3SCGAMA2) «0V2ZE*SIN (GAHA2) 
4 * S I N (G AHA1 1 

nv?7  = OV2XE*SrN(r,AMAi)  *OV2ZE*OOSfGAM4t> 

3V’E  = SQPT(3VZX'*.*2»0V2YE*»?*1V2ZE**2I 
■»V2XY?  = SaRT  <3V’X**2»0V2Y**2OV2Z**2> 

C 

C M34  COMPUTE  THE  THRUST  DIRECTION  AN5.ES  FOR  240  BURYI 


M. 
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®HIA=ASIN(0¥?7/OY2) 
®HI3*ASINO¥?YMO¥2*COS(PMIM)  ) 
T*0¥2K,Lr.  0. 1 9HI3*-(PI*PHIJ> 


NON  PRINT  RESULTS* 


PRINT*," 
PRINT*," 
PRINT*," 
PRINT* 
PRINT* 
PRINT*," 
»**  ",SY2 
^nr* 
print* 

PRINT*," 
PRINT* 
PRINT* 
PRI NT*  »" 
PRINT* 
®RINT*,“ 
PRINT*," 
PRINT*  ," 
PRINT*," 
“RINT* 
PRINT 
PRINT*," 
PRINT*  »“ 
PRINT*  »" 
**^I  *¥T* 
^IIT* 
®RI NT* »" 
PRINT*," 
®RINT*,“ 
PRINT*," 

print* 

PRINT*," 

®rtnt* 

PRINT*," 

®RINT* 

PRINT*," 

PRINT*," 

•»a» I MT * ,“ 
PRINT*," 
PRINT*," 
PRINT*," 

print* 

"RINT*,“ 

PRINT* 

PRINT*," 

PRINT* 

#U1T* 

PRINT*," 

PRINT* 

“RINT* 

PRINT* 

PRI NT* 

PRINT* 

PRINT*, “ 


‘XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 

•THE  TARGETING  RESULTS  IN  THE  LOCAL FRAME  ARE*  * 

‘XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 


• F3R  X PLANE  CHANGE  ANGLE  CSf ?)  DURING  THE  rIRST  BURN  OF 
?*R3,"  **  OEGREES,  THE  RESULTS  ARE*  " 


"XXXXXXXXXX  FOR  IH»ULSJ¥£  TRANSFER  XXXXXXXXXXXXXXXXXXXX" 


"TNE  TRANSFER  ORBIT  CHARACTERISTICS  ARE*  “ 

"THE  ECCENTRICITY  0r  THE  TRANSFER  ORBIT  IS* 
“THE  TRANSFER  ANGLE  IS*  ",CA*R0,"  OEG.  " 
“Y1T * ",¥1T,"  FT/SEC.  “ 

"¥2T«  ",¥2T,“  FT/SEC.  - 

IT*. *T HE  R2T  COHPUTEO  ANALYTICALLY  IS*  * 
“R2TX*  ",  R2X/CNMF,"  N, M.  " 

“R2TY*  ",  R2V/CNMF,"  N. N.  - 
"R2TZ«  ",  R22/CNNF.-  N.  H.  “ 


THE  COHPUTEO  ANGLES  ARE*  - 
THE  ANGLE  SY1  IS*  ",SYl*R0,"  OE5. 
THE  ANGLE  SYS  IS*  “,SY3*R0,"  DEG. 
THE  ANGLE  SV4  IS*  “,SYH*R0,*  OEG. 


“THE  TARGETING  INFORNAT ION  IS*  " 

“THE  rI»ST  BURN  INERTtAL  THRUST  DIRECTION  ANGLES  ARE* 


PH  1 1 * 
PHI?* 


'»  PHI1*R0,  1 
•,PHI2*R0,' 


"THE  2 NO  BURN  INERTIAL  THRUST  DIRECTION  ANGi.ES  ARE* 


PHIS* 

PHIA* 


•,PHI3*R0,' 
*»  PHIA*R9»  ‘ 


" THE  T3F*  IS*  “,T0F« 


xxxxxxxxxxxxxxxxxxxxxxxxyxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx" 


1 
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‘ I 

• I 


- 


®RINT*,-TNE  FDL LOWING  RESULTS  ARE  TO*  — FINITE  THRUST—  SURNSt  " 
PRINT*, «••••••••••••••••••••••••••••••••••*•**•••••••••••••••••- 

PRINT* 

PRINT* 

PRINT* 

269  RETURN 
END 

SU3R0UTINE  ISUESS 
C 

C XXXXXXXX XX XXXXXXXXXXX XX XX X<X<KXXXXX<K<X<XK<XX<XXXXXXXXX<<KXXXXX XXX 

C THIS  SU®RO JTI Nr  GENERATES  THE  CONPATI3LE  TRANSFER  TRAJECTORY  FOR 

C TNI  rIXE0  VE.OSITV  CHANGE  INCREMENTS  9T  A RES JL A FALSI  ITERATION 

C or  0X2 (SY1) • SECONO  BURN  DELTA  VELOCITY  A3  A rUNCTIDN  OF  THE  FIRST 

C DJRN  FLIGHT  »ATH  ANGLE,  UNTIL  ITERATED  DV2  HATCHES  ACTUAL  0V2. 

C XXXXXXXXXX<XXXXXXXXXXXXXXX<XXXXXXXXXXXXXXX<XXXXXXXXXXXXXXXXXXXXXXX 

• c 

CDHHON/TT/EHJ,®I,CNHF,RO,l)R»RlTfVlTfRtS,Vi:fR2T,V2Tf R2C,V2C,0V1, 
SOVJ.TNUl, TNU> ,i,CA,AT,SYl,SY2,SY3,SfL,I?N, THETA 
REAL  I2N 


C 

C 

C 


1 


6 

in 


c 

c 

c 


c 

c 


• I 


COMPUTE  HAXIHJH  POSSIBLE  SY1 

SYlHAX*ATA<MDVirVlS> 

PRINT* 

“RINT* 

ICD*0 
CRI T *• 01 
SY1«-,2*0R 
D5Y1*. 2*0R 
Sft«SYl*OSTl 

Tr(SYl.GE.SYlHSX)GO  TO  20 
GD  TD  10 

SYtN«SYlC-CCSTlC-SYlPl/(0V2C-0V2“n*(DV?C-DV2» 

SY1»SY1N 

91*ACOS(CDS($Y?)*COS(SY1>> 

OJANTO*(Vi:/DVn*SIN<81) 

Ie { A9S I QUANTO)  ST.l.tGO  TO  2 

A* ASIN  <QUANTf  I 

9*»I-A-B1 

0JANTt«0Vl»*>fVlC**2-2.*Vi:*DVl*C0S(3l 

TF(DUANT1..T,0  i GO  TO  2 

VI T* SORT (OJ ANTI  1 

Er*( VlT**2l /?• • EHU/R1T 

DJANT2*2.ME1U'R2T*ET> 

Tr ( 0UANT2* .T • 0 ) GO  TO  2 
V»T*5Q»T(DJANT?» 

AT«-'HU/(2.*ET) 

H«R1T*V1T*;0S(SY1» 

®*H**2/EHU 

E«SDRT<1.-VAT> 

CHECH  TO  ASS  JR*  VALIO  TRANSrER  0R9IT* 

Tr ( ®/ ( 1 • *E> • CT  R1T)G0  TO  2 
ir<P/Cl.-EI  ».T.  R2T)GO  TO  .2 
OJANTM«{ V2T*»2fV2C**2-OV2**2t /(2  *V2T*V2C> 
ir<  A9S  (O'JANTWI . GT.1JGO  TO  2 

CONPJTE  TRJE  A'lOMOLIES  AND  CENTRAL  ANGLE! 


. M.. 
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c 

TNJ1«AC0S((P/R1T-1.)/E> 

TNJ2*AC0S(<P'R?T-1.)SE) 

CA*TNU2-TNJt 

C 

C CHECK  FOR  »IANAR  OR  NON-COPLANAR  TRANSFER! 

C 

Ir( THETA. £1.0.) GO  TO  95 

C ND4  COMPUTE  TH*  RESULTANT  PLANE  CHANGE  RE3JIRS0  AT  SECOND  BURN! 

C 

6J1NT3*(SIH(3Y’>  /SINITHETAD’SINOAl 
IF| A9S (QUANTS)  • GT.l.) GO  TO  2 
Q«ASIN(QUANTSI 

0«2.*ATAN<(C3S<.5*<CA-0>nMC3S(.5*ICAO)> *TAN< .5* (TM£TA*SV2> >>> 

SYi»«PI-C 

GO  TO  56 

55  5YA»SY2 

56  r)'JANT4*H7CR2T*V2T» 

TF< A*S (OUANT4)  GT.l. ) GO  TO  2 
SY3*ACOSOJANT.» 

IF(SY3.LT.(20.»OR>>OSYl».ai*3R 

C 

C N34  COMPUTE  0Y2  IN  THE  E FRANEt 

C 

V2C*«0. 

Y?CYsV2C*03S<SY4> 

V2CIsV2C*StN<SY4> 

Y2r««Y2T*SIMCSYS» 

V?TY«Y2T*C3S<SY3) 

V2T7*0. 

OY2K-»Y2C«-¥2T< 

’W2YE»V2CY-V2TV 
9Y77£»  V2CT* Y2T? 

■W2S*SQRT<  JV2XE  **2*0V2VE**2OV22E**21 

Is (ICO. GT. 0 ) 53  TO  5 

TFOY2E.6T.OY2I  50  TO  7 

«?Y1P«SY1-.**3R 

SY1C«SY1 

OY2COY2E 

IC3*IC0*1 

r,3  T3  6 

7 9Y?»«0V2f 

60  TO  1 

5 9Y2N-9V2E 

IC(A9S(OV2NOY?) . LT.CRIT) GD  TO  11 
*?Y1»«SY1C 
SY1C-SY1N 
3Y?P*0¥2C 
3Y2C-0V2N 
63  T3  6 
2 60  TO  1 

20  PRINT*,"F3R  SY»«  ",SY2*R3,"  NO  IH»JlSIYE  TRAJECTORY  IS  POSSIBLE" 

6 PRINT* 

CA«0. 

63  T3  169 
11  CONTINUE 

PRINT* 

PRINT* 

3Y2*0Y2E 
169  RETURN 
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SUBROUTINE  CALEUN(N,r,P) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

tnis  routine  contains  the  nonlinear  equations  ihit  conditions) 

to  9*  SOLVEO  9*  NS01A. 

XXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

external  s.ope 

OINENSION  f <N),F (N) ,X (6),  2(5) , SNS (&,3) 

REAL  H0.MO,H01,NOl,N02,HO2,icl,HF2,I2, I?N,I1A,II 
CONNON/TVEHJ.PI.CNNF.ROOR,*!!,  W,RtC,N13,R?T,V2T,R2C,*2C,D¥l, 
*•>¥?,  TNU1,  TNU2,  ",CA,  AT,  SY1,SY2,  SYS,  STL,  I2N,  THETA 
CONNON/SLC/2 

COiNON/SLM01,N0  2,Mri,HF2,Tl,T2,N01,N02,T9A,T3ll 
CONNON/TE/X,SNS 

COHHON/OOAO/  R?CC,Y2CC,R2CA,R2CP,I2,E2CS,IIA,I1,TA,026 

C3NN0N/CH/NAXPL 

T.O. 

SF»l. 

ST»i.  0E5 

»(2)»Y(2)*SF 

2(1I»Y«3)»SF 

?u)»yo»)*sp 

2(J)«Y(5)*ST 

XI1)«R1C 

X(?)»l. 

X(S)M, 

XIV) *0# 

X(S)M1C 

X(S)«9. 

TT>TNA/9tt. 

SAUL  SET (S, T,«,OT, SLOPE, 0,.T., 0,3) 

00  10  K*l«$l? 

CA.L  STEP(i, T,X,OT, SLOPE, 0..T. ,0,0) 

continue 

CONNOTE  INCLINATION  OP  TRANSFER  ORUTl 

MIX«X<2)*X(S)-<<3)»X(9) 

MlY«X(3)*X(i»)-X(l)*X(*| 

N12»KC1)»«C9)-X<2I*X(RI 

Ml«S0RT(HlX*»2*mr**24HtI*«2) 

I1A*AC0S(N12/N1) 

3rs(?(9)-(T9A»*BO))/102M. 

CA.L  SET (S,T,X,OT, SLOPE, 0,.T., 0,0) 

00  20  «*1,19?<»0 

CA.L  STfcP(i,T,»,OT,SL0PE,0,  ,T. ,0,0) 

CONTINUE 

•)T«T09/912. 

ca.l  set (s,t,<,ot, slope, o,,t,, 0,0) 

00  30  K»l,*12 

CA.L  STEP(S,T»*» OT, SLOPE, 0, , T, ,0,0) 

CONTINUE 

CONP'JTE  RESULTING  NISSION  0R3IT  ELENSiTS* 

R?CC«SORT(<(l),*2»X(2)**2»X(3)**2) 
x?CC*SORT  (<  ( '■»)•*  2*X  (5)**2M(S)**2) 


140 


o 


t 


M2X«XC2»*X<6>-XIS>*X<5> 

H2V«X<3»*XU>-X(1>*XC6> 

H2I»XCH*X<5»-«<2»*X<6» 

*42»S'm»H2<*»2»H2Y**2*M2Z**2l 

XM«S<)RT(H2X**2*M2V**2> 

ir(rHer«.£).o.)6o  ro  uo 

o*;«aco$<-Y2Y/<H> 

IPIH2X.CT.0.>2?G«2.*PI-O2G 
kd  20YTINUE 

I2*020S«H2Z/H2» 

?9*sX(l)«X<4t  »x<2)*xm*xm*x<6> 

9<T«Y2CC**?-EYJ/R2QC 

SXMt.  /EHU>M9<T*X<l>-ROY*<<VI> 

EY«I1./EHJ)MB<T*XI2>-RW*X<5I  » 
E?«<1./EMJ)M9<T*X<3»-R9Y*<<5»> 

?2:C*SQRT(EX**?»EY**24E2**2» 

»2M2**2/£YU 

r2:p«p2/(i.»E2:ci 

r2:i»p2/u,-E2::> 

OUX>XU)-U2 
ORIY-XIEI 
^tr«xc3>  • 

■m«SQRTOU<*»20RlY**20m**2» 

TM»C0S<<R>C:**2*R1C**2-0R1**2>/<2.*222C*«C>) 

c 

5 Pnsr  COHP.IT-  » UNIT  VECTOR  tY  R2.H2,  AYJ  Y2C  OIRsCTIOHS 
C 

UY2X«X  ( D/Z223 

'«2Y»XI2)/?25: 

•JR2Z«X<3»/«22 
0H2X*H?X/H’ 
tM?Y*H2Y/HZ 
IM2?« H2Z/H2 

tlY?X«UM2Y*JR?Z-UH22*UR2Y 
•JY?Y*0H2Z*  JR*X-'JH2X*0R22 
'IW?2*UH2X*JR2Y-UH2Y*UR2X 
C 

2 N3Y  20HPUTI  Y3YIN6L  Y2  C3RRESPON3INS  TO  CURRiYT  R2  3IRE2TI0H. 
C 

Y2YX«Y2C»'J/2X 
Y2Y  Y*Y2C*'JY2Y 
Y?Y7»V2C*t>/2Z 
C 

C COY»UTE  P<Y»  -3*  CHOICE  3P  HIT  CONDITIONS* 

C 

T*<*«XPL.ST.0»;O  TO  50 

c 

<Mi>«<IlA-CIt*SV2)»Ml.EJ» 

‘■(2»«(B2c:-r?:»/(i.e5» 
p<  5) « Y2CC-Y23 
f<4>«E2CC*(1.E2I 
P( 5) » I T2N-I2) * ( I «E5) 
r.3  TO  60 
C 

90  CONTINUE 

«,d)»Y2MX-<(U  . 
rT2>»Y2NY«X (3) 
e<  J>  »Y2NZ“< ( St 
*<.>«EPCCM1.E?» 
c( ») * ( I2H*t  21  * ( 1 • £5) 


i 
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0 


c 

60  tET'JtN 

*19 

SU9S9UTI1E  S.9PE(N»T,X,0X) 

C 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C THIS  SOUTH*  391TAI1S  THE  EQJATI31S  0*  iOTIOM  TO  BE 

C IlTESSATEO  8T  SET/STEP. 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

9HE1SION  X 1 1>  , OX  ID  » Z (5> 

SEAL  10,10,101,H01,H0  2,H92,H*ltHF2,I2tt?»t 

C91131/SL/10l» i02,1Fl#MF2,Tl,r2,191»132# T§»*TIB 

C91101/SLCFZ 

<riJ«l.fc07&»*'16 

TF.T1  t TS-T2 

1*1  T.  ST.  T9AI  Tr*0« 

Tx*re*cos<zc?))*coscz(in 

TT*  TF*COS(Z  (Ml  *SI1  (Z  (II) 

TZ*TF*SHIt(2ll 

10*101 

10*101 

TH£*T 

IFlT.lT.  <Zm-TBB))GO  TO  2 
TX«TS*CO$<m>»  •COSCZ<3>> 

Tr«TS*COS<Z<V))*SINCZ<3» 

TZ«TS*SIN<IU>» 

10*102 

10*192 

TI1s*T-IZ«5»-T3B) 

2 301TI1UE 

S*53ST  (X  (1) **2*X  (2> *• 2+XI3)  **21 
9X1  1)  *X(4) 

OX(2)*X<5> 

9X(3>*X<6> 

OX  I '«>  «“E10*X(  1)  FS**3*TX/I  19 « ID* TI1EI 
OX(5>*-E1J*X(?»  FS**  3* TV/I  MO-19* THE) 
9XI6)»-EH.I*XITI/S**3*TZ/C10-19*TI1E> 

SETusi 

013 

SJSS3UTI1E  ETSOSSIPMIl.PMIZ.siIStPIlVfTO*) 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
*ms  SOUTHS  3OMP0TES  INSERTION  £RSO?  SEMSITI/ITIES  OUE  TO 
T1XJST  XESroT  iISAlIQNHEIT. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


"XTESIAt  S.OP- 

OH- IS  ION  <1rl  5*  »SlSI6t3>  ,XIS>  »Zt5) 

SEAL  10tM9t101» 101»  NO  2»M32» 1rl »1c2» 12 t 121 

001191 /TT/iU,»I  ,CM1»»X0»'9S,TlT,nT,?i:tVi:,?2T,W2T,P2C,X2Ct0Vi, 
5'W2tTiui,Hj?tr,CA,AT,SYt,sr2,srT,sv<»,r?N,risrA 
39ini/SL/101*l0  2fiFt  »NF2»  Tl»T2»19i t132»T9A»T3B 
09H01/TF/X  « SIS 
991101/SL3/Z 
1*0 

TS*,00i 

ST9S-:  THE  TATS-TEO  MOMXMAL  S2C  AID  V23» 


no  10  J«i.S 
XSF ( J) «X(JI 


I i 

I ‘ : 


t 


coib'jte  unit  vectors  in  the  iser^ial  using  thrust  angles* 

'J1XOOSCPSI2*  *;0SCPHI1» 

U1V«C0S<PMI2»*SIN<PHI1I 

U1Z«3INCPH12» 

U2X*C0SCPHIAJ  *OOSCPHIS) 

U2V«COS(PHt Al  •SINTPHI SI 
U2?«SIN(PHIA> 

COOROINATITE  ME  UNIT  VECTORS  IN  THE  SISALIGNED  PROSE 
FOR  ALIGNSENT  ERROR  ABOUT  X AXIS* 

'nxs»uix 

inVH*UtV*33S<TRJ-UtZ*SIN(TS» 

Ut?1«Uir*SIN(TR> ♦UlZ*COS<TR> 

'l>XN«U2X 

'J?VH«U2V*30SCTR>-U2Z*SIN(TR> 

U2ZH*U2Y*SIN<  TR) ♦U2?*C0S<TR1 
GO  TO  20  • 

*3R  ALIGNSENT  ERROR  ABOUT  THE  V AXIS* 

tJtXN«UlX*33S(TR>  Z*SIN(TR) 

U1VS«U1Y 

'J1ZN«-U1X*SIN<TR»«-UIZ*C0SCTR> 
tJ2XS*U2X*33S(  TR)  ♦U2Z*SINI TR) 

’)?VS»U2Y 

'J?ZN»-U2X*3IN<TRHU2Z*COS<TR> 

60  TO  20 

FOR  ALIGNSENT  ERROR  ABOUT  Z AXIS 

'J1XS«U1X*33SCT>)-U1Y*SINCTR> 

UlXS«UlX*SrMirRMUlT*COSCTR> 

IJlfH-UtZ 

•l?XH«U2X*33S<m-U2Y*SIN<TR) 

'/Zf  S«U2X*SmTR>*U2Y*C0S<TR> 

•I2TS«U2Z 

COMPUTE  THE  SIS ALIGNEO  THRJST  ANGLES* 

»HI2M«ASIM<UtZS» 

*»■*!  lHsACOSIUlXS/COS  (PHI 2H)  > 

»HIWM>AS1S<U’ZX) 

•*HI  3H» ASlNf  UJVH/COS  CPHIAH)  ) 

I*<  JZXH.LT.O.  I »HI3H»-«PI*PSIJN> 

T<  l ) *PHI IS 

”»I?)«PHI2S 

»m»PMISN 

’UMPHIM 

TOI-TOF 

GO  TO  IBS 

COS»UTE  RZC  AS0V2C  USING  SISALIGNED  THRUSTS* 

TO. 

xm»Ric 


143 


X<?|*0. 

xm«o. 
xm*o. 
xm  «vic 

X(i)»8# 

9T«T9A/S12. 

CO.L  SET(S,T,X,DT,SLOPE»0,.T.,9,0> 

•>3  70  K«l,512 

COuL  STEPtS  »T tt, 9T»  SLOPE* 0*  . T,  , 0*  0) 

70  CONTINUE 

OT*«TOF  -<T94»T99>>  /If  EM). 

CO.L  SET(5tT,*,3T,SLOPE,Of.T.,D,3) 

30  SO  X* I* 10200 

CO..L  STEP<itTfX,0T,SL0PE,3,.T.,0,3> 

80  CONTINUE 

3T*T09/512* 

0«LI  SET  <*s,T,<,OT, SLOPE ,0,.T.,  9,9) 

03  90  K«l,512 

CO..  STEP(S»T»X»OT» SLOPE* 0« • T, ,0*0) 

90  CONTINUE 

C 

C COMPUTE  SENSITIVITY  COEFFICIENTS  BY  C3LUHNI 

C 

I*I*1 

SNS<l,I>*C<<l»-XMF(l>  >/607S. 115*90 
SN5C2*  I>  «€<€»>  -XNF<  2)  > /607S. 115*38 
SN3<S,I1«(<(5>-XNF(S) 1/6079.115*89 
SNS(*.I>«Xm-XNP(*> 

SNSIJ.I»*XC5>-«NFCS» 

*NS(S.I>*X16»-XNF(8> 

TrCI«  EQ« 1)  63  TO  2 
T*<I. 50.2)10  T3  3 
RETURN 
5X3 

SD9R3UTINS  9 JRN1  CPMI1  »PHI  ?»  PNI 3 *0F, TO  »RlC  » VlC»R2Ct  Y2C) 

C 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

r.  Txrs  sunpcjttn-  computes  tne  *irst  burn  the  (T»i>  *3*  C3»lona* 

C T»ONSFERS,  0N3  THEN  COMPUTES  THE  OCTJOl  23RRESP0N3ING  INERTIAL 

C TiTUST  OIRECTION  ANCLES  PR3N  T91  ONO  S»":iPIE3  OR9ITOL  ELENENTS. 

C XVXKXXXXXXRKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 

C3NN9N/TT91/OIl,AI2#01,02*.82fT22,TEST 

*-0L  tl.I2.LI2 

’I « 3. 1 •»159?5S3*890 

*3*191. /PI 

■>*•»!/ lit. 

NT*»9. 

crit*.  (id 

ITIR.f 

123*0 

Ti*ott 

T2*AI2 

I'lTIST.OT.O.) .0  TO  IN 
•»NI?*-PHI2 
•NI*»-PHI* 
lb  CONTINUE 

Ml* V1C/R1C 
M?» V2C/R2C 

*5R»l2.*Pt>/(MSIMl-H*>> 
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print* 

PRINT*, “THE  COPLANAR  SYNODIC  PERI00  IS»  ■*,» ER,"  SEC*  " 

PRINT* 

PRINT* 

U02«L02-O2 
T*0  • 

C 

C CONPJTE  Til 

C 

go  to  10 

5 TN*TC-(<TC-T»>/CTAC-TAP>>*fTA5-TAI 

T*TN 

io  tf*wof 

ri:p*ric*cos<ni*t> 

RtCa*RlC*SIN(Nl*T> 

R2CP*R2C*C0$<J02*H2*TFI 

R?Ca=R2C*SlN<J02*W2*TF) 

’TSlRP*R2CP-RtC» 

6S„Ri*R2ca-Rt:a 

RlC«SlRT(RlC»**2*RlMP*tt 

R2C«SQRT(R2C»*»?*R2CQ**» 

nE.R«SQRT(0£:.AP**2*DELRQ**2> 

TAA«ACOS<(Ai:*»2*R2C**2-DE.R**2>/<*.*AlC*A2CI) 

T*(I3O.GT.J)S0  TO  5 
TF((TAP.GE,T».lNO.TRA.LT.TA))SO  TO  0 
60  TO  7 
B TP*T-OT 

TC*T 
Tt:<Tiu 
130*1 
60  TO  6 

7 ta»«tra 

T*T*OT 
60  TO  10 
5 TAN«*RA 

If (ASS  (TAN»TA) . IT.CRITIGO  TO  15 
TP*T3 
T3*TN 
TA»*TAC 
ta;*tan 
SO  TO  6 
IT  -51«T 

T3?«T61*T2! 

YI*N1*TB1 

n 

C N04  COMPUTE  THE  THRUST  DIRECTION  UNIT  VECTORS  IN  THE  LOCAL  FRANC 

C 

xsi*cos«PHti)*:oscPHi2) 

YS1»SIH(PHI1)*:0STPMI2) 

’Si*SlN(®HI2) 

<S»«COS(PHIJl*?OS(PHI*l  • 

YS»«"?IMCPHI3)  *:0S(PHIA) 
t»?»sin<phiai 
C 

f NON  COOROINATI-'E  THESE  THRJST  VECTORS  IN  THE  SEOCENTRIC-EaUATORIAL 
C rRANE  FOR  CORRESPONDING  RlCTTOl). 

C 

<t*(OOS(O))*;O3(Vl>*SIN(Ol>*0OS(Il>*SrN(Vl))*XGl*(C0S(Ol>*SIN(Vl»- 
S'lNT  Oil  *605(11)  *C0S(V1)  )•  YilM  SIN(Ol)  • 5INU 1) ) *TGl 
Yl«(SIN(Oi» •:0?(Vl)-COSI31)*C0S(Il)*StM(Vl)  )*<G1*(SIN(01)*SIN(V1>* 
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6CDS(31)*C35(I1)  *COS  CN11)*YGIM«SINCI1»  *COS(01»  > *251 
Tl*(-SIN(I1)*5TN<  VI  >)*XGlMStN<It>*233(  VI  )l*YIl*  1305(111  1*761 

x2*(cos(on*:osm>  ♦siN(on*:os(it>*siN(vii)*xG2*<:3scoii*siNm»- 
SSIN(01>*C35(lt»  *C0S(V1>  »*YG2*<SIN<01>*3IM<I1>>*7G2 

¥2*<sin<oi»  *335(vi>-cosm>  *:3nn>  *sinivi»  >*<g2*(sincoi>*sin(vi>* 
SCOSOi ) *035111)  *C0S(V1>>*  VSZM-SINdD’CnS  com  *ZS2 
?2*(-SINUt>*SrN<Vl»»*XG2*(SIN<It>*33S(Vin  •Y32*(303(Iin*ZGZ 

c 

C NO*  COMPUTE  THr  ACTUAL  INERTIAL  THRUST  AN5.ES* 

c 

»HI2A*ASIN<Zt> 

“HIIA=ASIN(Y1T30S(PHI2A>) 
in  (X1.AN0.Y1I  LT.O.)PHIlA«-<»I*PNItA) 
incXl.LT.l.I.SNO.  (YI.CT.0.>>PHI1A»<»I72.  *>HIU) 

®HI 4A»  ASIN(  Z21 

bhI3A=ASIN<Y2/:3S<PHI4A>> 

in  (X2.AN3.  T?>.  LT.0.>PHI3A«-(PI*PHI3A) 

in  (X2.LT.3.J.1NO. (Y2.GT. 0. >>PHI3A*(»I/2.*»HI5A> 

PRI NT* 

PRINT* 

ITER*ITER*l 
PRINT*, “ •****•  * 

PRINT* 

PRINT*, “P3R  NOMINAL  MISSION  START  TINE  NO.  ",tTER,"  TARGETING  IS*" 
PRINT* 

PRINT*, "1ST  3JRN  POSITION  ANILE  * “,V1*R0,"  DEG,  • 

PRINT* 

PRINT*, -T31", ITER, “«  “,TB1,"  SEC.  " 

PRINT*, “T32“,IT£R,"«  “,T32,“  SES.  * 

PRINT* 

PRINT*, “THE  THRUST  DIRECTION  ANGLES  ARE*  " 


PRINT*, “PHtl»* 

",PHI1A*R0,“ 

OEG. 

m 

PRINT*, “PHt2*« 

“,PHI2A*R0," 

DEG. 

m 

PRINT*, “PHI3A* 

",PHI3A*RD,“ 

0E3. 

m 

PRINT*, “PHI41* 

",PHI4A*R3,“ 

OEG. 

m 

T«T*PER 

Is ( ITER. LT • 5>  50 

TO  15 

RSTJRN 

ENO 

SUBROUTINE  R'JBNN  (PMI1, PHI  2,  PHI  3,  PHI  A,  3 23,  R1C, V1C«R23, V2C) 

C 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C THIS  SUHP3JTIN:  COMPUTES  THE  rIRST  3URN  TINE  «N0  CORRESPONDING 
S THRUST  OIRECT I3N  ANGLES  FROH  TB1  ANO  SPECIFIES  ORBITAL  ELEHEMTS, 

C C3R  N0N-C3»L»N»R  TRANSFERS. 

C XXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXX 

c 

C3HN0N/TTB1/AI1 ,AI2«01,02,L02,T22,TEST 

REAL  I 1,I2,L02 

PI«3.14t53»5S3SM$ 

ROMNO./PI 

nR*®I/l»0, 

It* All 
I?«AI2 
TTER«0 

ip«t:st.gt,o.)go  to  i% 

BHI2.-PHI 2 
PHIA*-PMIA 
14  CONTINUE 

HI*  V1C/R13 
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M2«V2C/R2C 

COMPUTE  TSl 

111*02-026 
IF(Ul.GE.O.)  arc*ui 
tfiui.lt. o.»arc*(2.*pi-ui» 

TB1-ARC/W1 

“ER«(2.*PI)/Nl 

T9?*T81*T22 

MOM  COMPUTE  THE  THRUST  DIRECTION  UNIT  SECTORS  IN  THE  LOCAL  FRAME 

x3t«cos(PHtii*cos(PHi2> 

Y31 *31 N(PHI II  *C0S(PMI2) 

Y31*SlN(PNI2l 

XS2*COS(PHI3»*:OSCPMI%> 

YS?*3IN(PHr J»  *COS(PMIA> 

?S2*SIN(PHI4I 

N3N  C00R0INATI7E  THESE  THRJST  VECTORS  IN  THE  3E0CENTRIC-E0UAT0RIAL 
eRANE  FOR  C0RRES«»0N0IN6  R1C(TS1>. 


XI*  (COS  (Oil  *C3?  (VII  ♦SIN(31I*C3S(I1>*SIN(VI 
5SIH(01)*COS(Il)*C0S(Vl) )*Y31*(SIN(01)*5TN( 

ri*(5iN(oi»*C33(vn-cosoi»*cos(ii)*5ry(vi 
&C3S (01 I *COS( It  I *COS (Vl> >*Y31*(-SIN(Il) *C3S 
Tl*(-SIN(Il)*STN(Vtl)  *XG1  *(SIN  (ID  • COS (VI) 
X2« (COS(01)*C3S(V1>*SIN(31I*COS(I1)*SIN(V1 
5 SIN (31 >*G35 (III *COS (Vll >*Y32M  SIN(01>  *SIN( 
Y2*(SIN(0H *33S(V1>-C0S(31I*33S(I1)*SIN(V1 
3 COS (31  ) *C33 (I1)*C0S(V1) )*VS2M-SIN(Il)*C05 
Y2*(«SIN(III*SIN(V1)I *XG2*(SIN(I1) *C35 ( VII 


)»*<Gl*(C3S(01l*SINm>- 
1 1)  I *ZG1 

» >*<G1*(5IN(01»*SIN(V1>* 
(Oil  I *731 

l*Y3l*(C05(Ill)*  261 
I)*XG2*(C3S(01)*SIN(V1>- 
1 1>  > *ZG2 

>>*<G2*(SIN(01I*SIN(V11* 

(01I»*7G2 

)*Y32*(C03(I1)I *Z62 


N3N  COMPUTE  TN"  ACTUAL  INERTIAL  THRUST  AN6.ESI 

*HI2*»ASIN( Zl) 

PHria*ASIN(riFC3S(PHI2A)> 

IF((X1.AN3.  Til  LT.0.>PHIH«-(PI*°HI1AI 
Tc((X1.LT.O.).ANO. (Vl.GT. 0. > >PHIlA*(»I/2.  »“HIIA) 
“HH»A«ASIN(Z2> 

•>MISA*ASrN»  Y>FCOS(PHIfcA>) 

IF(  (X2.LT.3.I,ANO.(Y2.6T.O.  >I“HIU«-(PI*“HI3A) 

TF((X2.AN3.Y2)  LT. 0 . > PHI3 A*-( »I*»HI3A> 

Ic( (X2.LT.0.I .MO. (Y2.GT. 0. II PHI3A* (“I /2. *»HI3A) 

“RINT* 

“RINT* 

“RINT*,-  • ••*•  “ 

PRINT* 

PRrNr*,-FO(  NOMTNAL  MISSION  START  TINE  TSl,  THE  TAR3ETIN6  ISI 
“RINT* 

“RINT* ."1ST  3JRH  POSITION  AN3.E  * ".Jt*R3,"  OEG.  * 

“RINT*, "THE  *• RIOD*  “,PER,“  SEC.  " 

“RINT* 

“RINT* 

ITER«ITER*l 

“RINT* ,"T9t“»ITFR,"«  ",T81,"  SEC.  " 

“RINT*, ITER, "»  ",T32,"  SEC.  " 

“RINT* 

TRI»TR1*PER 


14? 


»■ 


IF(ITER.LT.5>S0 

PRINT* 

PRINT* 

TO  10 

PRINT*,“TNE  TNRUST  OIRECTIOM 

ANGLES 

ARE* 

“RINT*  ,“PHI1A* 

~,PHI1A*R0,- 

DEG. 

m 

®RI NT*  ,“PNI2A* 

“,PHI2A*R0,“ 

DEG. 

m 

PRINT* ,“PNI3A« 

“,PMI3A*R0,- 

DEG. 

• 

PRINT*, -PHIR1* 

",PHIRA*R3,** 

DEG. 

• 

PRINT* 
p^r mt* 
re torn 

ENO 

Sl)3  ROUTINE  HRJSTF<RlCfV13,T0F,Cl,S2» 

c 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 'HIS  ROUTINE  OONPUTES  INSERTION  ERROR  SHNST TIi/I TIES  DUE 

C TO  thrust  nasnituoe  DEVIATIONS. 

r XXX XXXXXXXXXXXVX XXXXXXXXXXXX<XXXXXXXXXXXXX<XXXXXXXX<XXXXXXXX 

c 

EXTERNAL  S.OPE 

OUENSION  <N<S)  t X (6 ) , SNS(6,  3)  , Z (5) 

53iN0N/TE/X,SNS 

03NNON/SL/N01,102fMFl,HF2#Tl,T2,101,N02,T81,T38 

CONNON/SLC/Z 

REAL  N01.N02,  H-l,  HF2,M01tN02,NDlN,N02N 
CNNFs6076.U54S6 
IS<IP*0 
NDIN«HD1 
•O2N«N02 
T8AN*T8A 
TR3N*TB8 
T1N*T1 
T2N«T2 
30  10  J»1,S 
10  XN<J)*X<J> 

ti«ti*«.oi»*ti 

T?*T2»  x. oi»  *rz 
N01*T1/C1 
•02*72/02 
T8A*<N01-NF1|/101 
T33*<N02“Hr2l /ID  2 
100  T«J, 

X(1J«R1C 

X<2)*0. 

X(3»*0. 

Xf.»*0. 

X< J)*V1C 
X(i»*0. 

3T*TBA/512, 

CA.l  SET C5,T,X,0T, SLOPE, D,,T.,0,0l 
30  70  K-1,512 

SA.L  STEP(S,T,X,0T,SL0PE,0,.T.,0,D» 

70  oontinue 

3T«(T0F  -<T3*»TB3>>/102WI. 

3ALI  set d,T ,X,0-T, SLOPE, 0».T.» 0,01 
30  SO  K*1,10?AB 

OAlL  STE»(S,T,y,DT,SL0PE,D,.T.,0,0> 

00  CONTINUE 

3T«T98/512. 
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Cft.L  SET (5»T ,<»DT,SlOPE,Dt»  T. , 0»D> 

TO  99  K*1(512 

CA.l  STEPC>,TtX,OT,SUOPE,D,.T.  ,0,0) 

OONTTNUE 

iFtisxiP.Gr.eiso  to  z 

09**  (SORT  ( ( X ( 1 ) - XN(  1)  >**2*CXC2>-XNC2>»***MXm-XN(S>>**2>  >/CNHP 
OV*«SQRT  <<<<♦»*  KN(4>  I **2M<  (?) -XNC5II  »•*♦(<  (S)-XN<  ill  **E> 

T1*T1N 

T2*T?N 

T1*T1-  (• Oil *Tl 

T2*T2-  (• Oil *T2 

N01»T1/Ct 

ND2*T2/C2 

T34  * (M01-NF1  I /NOt 

TB3*CN02-HF2I/N0« 

IS<I®*1 
00  TO  101 
CONTINUE 

091*  (SORT ( (X  ( 1 ) - XN  ( 1)  >**2*(X<2>-XN<2) I CXCS) -XN (3) > **EI  I /CNNF 

OVN«SQRT(CXU>-XNf4))**2MX<i>-XNC5>l**2«-<<(6l-XN<6ll**2> 

•9INT* 

•9INT* 

•91 IT* 

•9INT*  .“THE  INSERTION  POSITION  AND  VE.DCTTV  E990R  SENSITIVITIES  " 
•91  NT*  ."F09  -TNRUST  NAG.  DEVIATIONS-  ,90TN  PLJS  ANO  NINOS  ARE*  " 
•91  NT* 

•91 NT* 

•9INT*  ,**09**  "♦  ORP»“  (N.N.I/CMX)  " 

•9INT*,"0V»«  -,OVP,"  (FT/SSC. > / {*1X1  " 


’♦ORP,"  (N.N.I/C+1X)  • 

■.OVP,**  (FT/SSC. >/I*tX> 
’» ORN,"  CN.N.I/C-IX)  - 

•,ovn,-  <ft/se:.)/(-ix> 


®9INT*,-0VN«  ",OVN,-  (FT/SEO 

•9INT* 

•9INT* 

•91  NT*  — — — — — — — — 

T1*T1N 

T2*TEN 

N0t*N01M 

N02«N02N 

TBA-T9AN 

T1J.T8PN 

9ETJ9N 

•NO 

SUBROUTINE  NTSpl A.N.N. AT) 

TIN* NS  ION  1 ( N.N) .AT CN.N) 

00  2 1*1 »N 
00  3 J*1»N 
1TCJ.I)*A(I,J> 

CONTINUE 
9ET J9M 
ENO 

SUBROUTINE  HNPfC A,3,C.H,KtN> 
OIN£NS ION  1CN.O  fO(<fN).C(N.N) 
TO  10  J*1,N 
00  10  I *1.1 
CCt'JI *0.1 
00  10  t«t,< 

SCI, J>  *CCI» J**i<I,L)*8CL,  J> 

RETURN 

ENJ 
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